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20.  ABSTRACT  (Continued). 

product  of  subsequent  model  development  and  experimental  studies,  Dissolved 
Oxygen  Subroutine  with  Anaerobic  Processes  (DOSAP) ,  is  presented  herein  and  its 
operation  is  described  in  detail.  The  DOSAP  model  is  presented  as  a  series  of 
seven  linked  geochemical  steps  (phases)  of  increasingly  anaerobic  (reduced) 
conditions.  Major  factors  influencing  the  rate  and  duration  of  each  of  the 
major  pathways  are  discussed.  The  DOSAP  model  has  subsequently  been  transformed 
into  a  mathematical  subroutine  for  the  CE  one-dimensional  reservoir  water 
quality  model,  CE-QUAL-R1. 

A  literature  search  revealed  a  paucity  of  information  on  various 
physicochemical  and  biological  rates  affecting  geochemical/water  quality  condi¬ 
tions  in  reservoirs,  although  abundant  data  .for  some  parameters  exist  for  other 
freshwater  systems  and  marine  environment s.^nA  series  of  investigations  con¬ 
ducted  in  250-t  laboratory  reaction  chambers jprovided  an  initial  set  of  physico¬ 
chemical  rate  data  for  DeGray  Reservoir  on  tie  Caddo  River  in  Arkansas.  Use  of 
large  reactor  systems  appears  to  be  a  feasible  approach  for  obtaining  the 
physicochemical  rate  data  required  to  use  DOBAP  on  a  site-specific  basis,  and 
this  method  is  being  used  to  develop  a  set  o ft  rate  coefficients  using  sediments 
from  CE  reservoirs  on  a  nationwide  basis. 

While  detailed  guidance  on  data  requirements,  sampling  procedures,  and 
rate  coefficient  selection  for  DOSAP  will  be  presented  in  a  later  Engineering 
Technical  Letter,  this  report  summarizes  procedures  currently  being  used 
to  obtain  data  for  the  subroutine. 
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PREFACE 

This  study  was  sponsored  by  the  Office,  Chief  of  Engineers  (OCE), 

U.  S.  Army,  as  a  part  of  the  Environmental  Water  Quality  and  Operational 
Studies  (EWQOS),  Work  Unit  3159^  (IB.2),  entitled  "Develop  and  Verify 
Descriptions  for  Reservoir  Chemical  Processes." 

The  work  was  conducted  during  the  period  September  19?8-June  1979 
by  the  Environmental  Laboratory,  U.  S.  Army  Engineer  Waterways  Experi¬ 
ment  Station  (WES),  Vicksburg,  Mississippi,  under  the  direction  of 
Dr.  J.  Harrison,  Chief  of  the  Environmental  Laboratory  (EL),  and  under 
the  general  supervision  of  Dr.  R.  L.  Eley,  Chief  of  the  Ecosystem 
Research  and  Simulation  Division  (ERSD),  and  Dr.  R.  M.  Engler,  Chief 
of  the  Ecological  Effects  and  Regulatory  Criteria  Group.  Program 
Manager  of  EWQOS  was  Dr.  J.  L.  Mahloch,  EL.  Some  of  the  foundation  work 
for  this  study  was  performed  as  part  of  a  Corps  of  Engineers  In-House 
Laboratory  Independent  Research  Program  (ILIR)  Work  Unit  111  06  entitled 
"Mechanisms  that  Regulate  the  Degree  of  Oxidation-Reduction  in 
Anaerobic  Sediments  and  Natural  Water  Systems." 

This  study  was  conducted  by  Dr.  D.  Gianni  son  and  Mr.  J.  M.  Brannon, 
ERSD.  Mr.  I.  Smith,  Jr.,  Ms.  P.  L.  Butler,  Mr.  G.  A.  Burton,  and 
Ms.-K.  M.  Preston,  ERSD,  assisted  with  the  laboratory  experimentation. 
This  report  was  written  by  Dr.  Gunnison  and  Mr.  Brannon  and  was  re¬ 
viewed  by  Dr.  K.  W.  Thornton,  and  Dr.  R.  L.  Chen,  ERSD,  and  Dr.  P.  G. 
Malone,  Environmental  Engineering  Division. 

Directors  of  the  WES  during  this  study  and  the  preparation  and 
publication  of  this  report  were  COL  John  L.  Cannon,  CE,  and  COL  Nelson  P. 
Conover,  CE.  Technical  Director  was  Mr.  F.  R.  Brown. 
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CHARACTERIZATION  OF  ANAEROBIC  CHEMICAL  PROCESSES  IN  RESERVOIRS 


PROBLEM  DESCRIPTION  AND  CONCEPTUAL  MODEL  FORMULATION 


PART  I :  INTRODUCTION 

Background 


1.  Many  Corps  of  E:  gineers  (CE)  reservoirs  experience  low  concen¬ 
trations  of  dissolved  oxygen  and  high  concentrations  of  products  of  an¬ 
aerobic  transformation  within  the  hypolimnion  and,  subsequently,  in 
waters  released  from  projects  having  bottom  withdrawal.  Adverse  impacts 
of  low  dissolved  oxygen  in  reservoirs  and  their  releases  have  included 
harmful  effects  on  aquatic  biota  such  as  fish  kills,  loss  of  project 
benefits,  high  operation  and  maintenance  costs  especially  as  a  result  of 
corrosion,  and  increased  costs  to  downstream  water  users  due  to  increased 
requirements  for  water  treatment. 

2.  Existing  modeling  approaches  do  not  adequately  predict  the 
many  biological,  chemical,  and  physical  ramifications  of  anaerobic  pro¬ 
cesses  that  occur  in  reservoirs,  although  several  models  of  pertinent 
individual  processes  do  exist:  DOSAG  of  the  Texas  Water  Development 
Board  (1970)  and  the  phosphorus  models  of  Jorgensen  et  al.  (1975)  and 
Kamp-Nielson  (1975,  1976).  Moreover,  present  techniques  give  an  in¬ 
effective  characterization  of  the  sources  and  extent  of  dissolved 
oxygen  depletion,  the  development  of  anaerobic  conditions  in  the 
hypolimnion,  the  mobilization  of  nutrients  and  contaminants,  and 
sediment  deposition  and  resuspension.  Managers  also  lack  the  capa¬ 
bility  to  evaluate  various  alternative  methods  for  correcting  such 
problems  or  to  determine  the  subsequent  effects  of  corrective  measures 
on  impoundment  water  quality. 


Purpose 

3.  This  study  was  undertaken  to  provide  the  information  for  the 
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development  of  an  anaerobic  model  which  could  then  be  used  as  a  tool  to 
deal  with  water  quality  problems.  This  report  presents  the  results  of 
an  investigation  on  processes  associated  with  low  dissolved  oxygen  and 
anaerobic  conditions  in  reservoirs.  This  review  has  led  to  the  formula¬ 
tion  of  a  conceptual  model  of  reservoir  anaerobic  processes  and  to  the 
development  of  a  descriptive  subroutine:  Dissolved  Oxygen  Subroutine 
with  Anaerobic  Processes  (DOSAP).  The  anaerobic  subroutine  has  sub¬ 
sequently  been  incorporated  into  the  CE  one-dimensional  reservoir  water 
quality  model  (CE-QUAL-Rl)  as  a  part  of  the  Environmental  and  Water 
Quality  Operational  Studies  (EWQOS)  Task  IC.  After  appropriate  testing 
and  evaluation,  the  subroutine  will  be  documented  in  subsequent  updates 
of  the  CE-QUAL-Rl  user's  manual.  This  report  is  confined  to  a  technical 
description  of  basic  changes  in  reservoir  water  quality  created  by  low 
dissolved  oxygen  and  anaerobic  conditions  and  to  the  conceptual  model 
that  was  formulated  to  describe  these  processes. 
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PART  II:  PROBLEMS  WITH  LOW  DISSOLVED  OXYGEN  AND 
ANAEROBIC  CONDITIONS  IN  CE  RESERVOIRS 


1*.  This  portion  of  the  report  discusses  problems  associated  with 
the  development  of  low  dissolved  oxygen  (DO)  and  anaerobic  conditions 
in  CE  reservoirs.  Consideration  will  be  given  to  depletion  of  dissolved 
oxygen,  to  anaerobic  nutrient  regeneration  and  contaminant  mobilization, 
and  to  the  current  strategies  used  by  CE  Field  Offices  to  address  these 
problems.  Keeley  et  al.  (1978)  discussed  these  considerations  in  detail 
Keeley  et  al.  (197&)  is  summarized  below  to  provide  a  background  for  the 
technical  discussion  presented  in  Part  III  of  this  report,  beginning  on 
page  13. 


Dissolved  Oxygen  and  Anaerobic  Conditions 


5.  The  impoundment  of  waters  often  results  in  the  development  of 
low  dissolved  oxygen  concentrations  within  the  impoundment  and  in  reser¬ 
voir  releases.  This  situation  may  result  in  the  accompanying  release  of 
undesirable  products  of  anaerobic  processes,  including  dissolved,  reduced 
forms  of  nitrogen;  iron  and  manganese;  reduced  organic  compounds;  and 
hydrogen  sulfide.  The  problem  is  particularly  relevant  because  of  the 
adverse  consequences  of  these  conditions  on  recreational  and  fishery 
resources  within  and  downstream  of  the  impoundment  and  because  of  poten¬ 
tial  increases  in  water  treatment  costs  incurred  by  users. 

6.  Many  reservoirs,  like  natural  lakes,  develop  thermal  stratifi¬ 
cation  (Figure  l).  Reservoirs  thus  affected  may  subsequently  develop  low 
dissolved  oxygen  concentrations  in  hypolimnetic  waters  because  the  strati¬ 
fication  serves  as  a  barrier  to  vertical  mixing  with  oxygenated  surface 
waters.  Natural  reaeration  of  hypolimnetic  waters  is  consequently 
minimal  until  destratification  occurs.  The  low  levels  of  dissolved  oxy¬ 
gen  and  the  anaerobic  conditions  that  may  ensue  are  the  result  of  a  com¬ 
plex,  integrated  series  of  physical  and  biogeochemical  processes  that 

use  oxygen-demanding  organic  matter  and  various  inorganic  substances. 

7.  Biogeochemical  processes  occurring  under  anaerobic  conditions 
differ  greatly  from  those  encountered  in  aerobic  environments.  Microbial 
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ATMOSPHERIC  REAERATiON 

- 

•  WARM  ISOTHERMIC 
EPlLlMNlON  •  ABUNDANT  OXYGEN 

•  WARMWATER  FISHERY 

•  WARM  TO  COLD  THERMAL  DISCONTINUITY 
METALIMNION  •  VARIABLE  OXYGEN 

•  Ml  XED  FiSHE  RY 


•  COLO  1SOTHERM1C 

•  OXYGEN  LOW  OR  ABSENT  INCREASED 
HYPOL'MNlON  CONCENTRATIONS  OF  SOLUBLE  FORMS 

of  contaminants  and  nutrients 

•  COLDVVATER  FISHERY  iF  OXYGEN  ADECuATF 

SEDIMENT  _  REGION  OF  MATER'AL 

A  DSORPT  i  ON  AND  RELEASE 


*  TYPICAL  VERTICAL  TEMPERATURE  AND  DO 
DISTRIBUTIONS  DURING  STRATIFICATION: 


Figure  1.  Thermally  stratified  reservoir  and  associated 
conditions  of  low  dissolved  oxygen 

activities  in  anaerobic  environments  yield  soluble,  often  highly  concen¬ 
trated,  forms  of  nitrogen,  iron,  manganese,  and  sulfides  that  can  produce 
toxicity,  taste,  odor,  and  staining  problems.  Moreover,  anaerobic  pro¬ 
cesses  also  release  nutrients  that  may  support  undesirable  algal  blooms. 
Finally,  extended  evolution  of  hydrogen  sulfide  in  reservoirs  can  result 
in  severe  problems  with  odors,  corrosion  of  outlet  structures,  and  toxic¬ 
ity  to  aquatic  organisms. 

8.  Several  approaches  have  been  used  to  deal  with  problems  associ¬ 
ated  with  low  DO  and  anaerobic  conditions  in  reservoirs  and  their  re¬ 
leases.  The  approaches  have  generally  been  very  expensive  and  have  often 
been  effective  only  for  small  water  bodies.  In  the  case  of  release 
waters,  existing  approaches  have  only  temporarily  or  partially  achieved 
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desired  water  quality  conditions.  Examples  include  injection  of  air  or 
oxygen  into  the  hypolimnion,  mechanical  destratification,  penstock  in¬ 
jection  of  oxygen  or  turbine  venting,  and  selective  withdrawal.  Figure  2 
illustrates  several  of  these  procedures.  In  most  cases,  inadequate  guid¬ 
ance  is  available  on  the  probable  severity  of  anticipated  problems  or  on 
when  to  initiate  and  when  to  terminate  corrective  procedures. 

9.  A  final  management  strategy  is  the  reduction  of  oxygen¬ 
demanding  materials  entering  reservoirs.  Reduction  may  be  accomplished 
in  several  ways:  diversion  of  waters  containing  oxygen-demanding  sub¬ 
stances  (Edmondson  1970);  appropriate  clearing  and  site  preparation  prior 
to  reservoir  filling;  or  filling  and  flushing  or  incrementally  filling 
the  reservoir  basin  during  initial  stage;  f  impoundment.  Alternatively, 
selective  withdrawal  procedures  may  be  used  in  reservoirs  to  avoid  the 
release  of  low  DO  waters  by  removal  from  various  levels  in  an  impound¬ 
ment  (Bohan  and  Grace  1973)  (Figure  2b). 


Nutrient  Regeneration 

10.  Excessive  nutrient  enrichment  of  reservoirs  can  promote  the 
development  of  nuisance  algal  blooms  and  the  excessive  growth  of  sub¬ 
merged  plants.  The  decomposition  of  these  plant  materials  may  produce 
fish  kills,  taste  and  odor  problems,  and  unaesthetic  conditions  (Keeley 
et  al.  1978).  Secondary  effects  that  result  from  the  presence  of  exces¬ 
sive  algal  growths  include  extreme  fluctuations  in  diurnal  oxygen  concen¬ 
trations  and  development  of  rough  fish  pojulations. 

11.  A  major  source  of  nutrient  input  into  many  reservoirs  is  run¬ 
off  from  agricultural  lands.  Excessive  inputs  of  nutrients  in  some 
areas  can  also  be  attributed  to  poor  management  of  forested  lands  and 
to  mining  activities  in  upstream  portions  of  the  watershed.  Nitrate- 
and  phosphate-rich  runoff  from  feedlots  and  cattle  grazing  near  the 
reservoir  are  other  sources  of  enrichment.  However,  even  if  nutrient 
inputs  are  reduced,  reservoir  sediments  under  anaerobic  conditions  may 
serve  as  sources  of  nutrients,  particularly  nitrogen  and  phosphorus, 
through  the  process  of  anaerobic  nutrient  cycling. 
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b.  Generalized  schematic  of  a  selective  withdrawal 
structure  for  reservoir  releases 


Figure  2.  Alternatives  for  problems  associated  with  low  DO  and 
anaerobic  conditions  in  reservoirs 


12.  Problems  with  excessive  nutrient  enrichment  and  DO  demands  in 
reservoirs  frequently  are  most  intense  during  the  first  few  years  of 
impoundment .  New  reservoirs  not  only  are  subjected  to  external  (alloch¬ 
thonous)  loadings,  but  also  contain  a  significant  internal  loading  and 
oxygen  demand  resulting  from  the  decomposition  of  inundated  herbaceous 
plants,  litter,  and  organic  rich  topsoil.  By  contrast,  the  decomposi¬ 
tion  of  woody  plant  material  such  as  trees  probably  is  not  a  significant 
cause  of  eutrophication  and  low  DO  because  of  the  very  slow  rate  of 
degradation.  Effective  preimpoundment  clearing  and  site  preparation 
policies  cannot  be  formulated  presently  because  of  insufficient  informa¬ 
tion  on  the  relative  importance  of  various  materials  as  sinks  for  DO 
and  sources  of  nutrients  following  initial  reservoir  filling  and 
inundation. 

13.  Many  of  the  watershed  activities  that  contribute  to  the 
eutrophication  of  CE  reservoirs  cannot  be  controlled;  however,  the  CE 
can  influence  the  effects  of  nutrient  regeneration  through  project  design 
and  operation  (Keeley  et  al.  1978).  Any  procedure  that  maintains  hypo- 
limnetic  DO  will  generally  minimize  the  impact  of  anaerobic  nutrient 
regeneration.  Restoration  techniques,  such  as  dredging,  that  remove 
nutrient-laden  sediments  from  the  reservoir  may  be  effective,  but  econom¬ 
ically  infeasible,  as  alternative  operational  strategies  that  prevent 
the  initial  stagnation  of  hypolimnetic  waters.  The  need  for  these  ex¬ 
pensive  corrective  procedures  might  often  be  avoided  if  adequate  tech¬ 
niques  for  predicting  the  nature  and  extent  of  anaerobic  nutrient  regen¬ 
eration  in  reservoirs  were  available. 

Mobilization  of  Contaminants 

14.  Contaminants  have  been  defined  as  certain  microorganisms  or 
chemical  substances  that  may  have  a  deleterious  effect  on  the  environ¬ 
mental  quality  of  a  given  aquatic  resource,  primarily  through  their 
adverse  effects  on  man  or  on  biological  food  chains  (Keeley  et  al.  1978). 
From  the  standpoint  of  anaerobic  conditions,  only  two  types  of  contami¬ 
nants  are  of  major  concern.  The  first  type  includes  those 
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oxygen-demanding  waste  substances,  particularly  organic  compounds,  that 
deplete  oxygen  from  water  through  autocatalytic  processes  or  that  are 
sufficiently  abundant  and  biodegradable  that  their  decomposition  will 
remove  oxygen  from  the  water.  The  second  class  of  contaminants  includes 
those  metals  whose  chemical  or  physical  forms  are  altered  under  anaerobic 
conditions  with  the  result  that  their  solubility  is  modified  and  their 
concentration  in  the  reservoir  is  increased. 

15.  Contaminants  are  a  source  of  problems  in  CE  reservoirs  pri¬ 
marily  because  certain  chemical  contaminants  such  as  heavy  metals  (e.g. 
mercury)  and  certain  pesticides  (e.g.  chlorinated  hydrocarbons)  are 
capable  of  undergoing  biotransformation,  and  often  bioaccumulation,  in 
food  chains.  Reservoirs  that  presently  do  not  have  immediate  contami¬ 
nant  problems  may  be  subjected  to  continual  low-level  contaminant  load¬ 
ings,  which  have  the  potential  for  being  accumulated  and  causing  signif¬ 
icant  problems  in  later  years  of  project  life.  This  is  particularly 
true  if  the  loadings  become  deposited  in  the  reservoir  in  locations  where 
major  biogeochemical  activities  may  alter  their  solubility  through  anaer¬ 
obic  processes.  Two  fairly  common  metals,  iron  and  manganese,  may  reach 
such  high  concentrations  in  anoxic  hypolimnetic  waters  that  reservoir 
discharges  become  discolored  and  cloudy  when  the  reoxidation  of  these 
materials  changes  their  properties  from  transparent  dissolved  to  colored 
colloidal  substances. 

16.  Three  general  approaches  are  currently  used  for  the  manage¬ 
ment  of  contaminants  mobilized  by  processes  active  under  anaerobic  condi¬ 
tions.  These  include  diversion  of  contaminant-laden  inflows,  removal  of 
contaminant-laden  sediments,  and  maintenance  of  aerobic  conditions  in 
the  reservoir  hypolimnion. 
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PART  III:  DEVELOPMENT  OF  LOW  DISSOLVED  OXYGEN 
AND  ANAEROBIC  CONDITIONS  IN  CE  RESERVOIRS 


17.  The  intensity  of  oxidative  and  reductive  processes  occurring 
in  anaerobic  sediments  and  natural  water  systems  directly  affects  CE 
Civil  Works  activities.  The  environmental  problems  associated  with  re¬ 
leases  of  anaerobic  waters  from  reservoir  hypolimnions  have  been  of 
special  concern.  The  basic  principles  of  DO  depletion  and  subsequent 
anaerobic  reactions  and  associated  environmental  controls  must  be  under¬ 
stood  before  a  predictive  capability  can  be  developed  for  anaerobic  prob¬ 
lems.  The  capacity  to  forecast  the  extent  and  intensity  of  anaerobic 
conditions  achieved  under  various  environmental  constraints  is  important 
since  the  intensity  of  anaerobic  conditions  will  determine  the  form  and 
mobility  of  many  chemical  constituents.  The  form  and  mobility  of  a  given 
chemical  constituent  will,  in  turn,  determine  the  potential  of  a  sediment 
or  water  system  to  exert  an  adverse  environmental  impact. 

Li terature  Review 


Reduction  process 
in  natural  systems 

18.  Natural  systems  are  considerably  different  from  rigidly  con¬ 
trolled  laboratory  systems  having  fixed  characteristics  (Clark  1972) 
since  natural  systems  are  constantly  in  a  state  of  flux.  In  addition  to 
changes  in  oxidation-reduction  potential  (Eh),  natural  systems  also  ex¬ 
hibit  constant  changes  in  chemical  composition,  hydrogen  ion  concentra¬ 
tion,  and  diversities  of  microbial  activities  (Hewitt  1933).  Oxidation- 
reduction  potentials  in  natural  systems,  however,  do  vary  between 
well-defined  limits.  The  upper  and  lower  limits,  respectively,  of  Eh 
in  natural  systems  are  determined  by  the  presence  or  absence  of  DO  and 
by  the  point  at  which  water  becomes  unstable  and  its  components  are 
oxidized  to  molecular  oxygen  and  reduced  to  molecular  hydrogen  (Garrels 
and  Christ  1956);  the  latter  point  is  approached  when  anaerobic  bacteria 
decompose  organic  matter  and  reduce  molecular  hydrogen.  The  oxidation 
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states  of  many  organic  and  inorganic  constituents  are  affected  by  the  Eh 
of  surrounding  water  and  sediment. 


Microbial  effects  and  bio- 
degradation  of  organic  matter 

19.  When  the  oxygen  supply  to  a  flooded  soil,  sediment,  or  aquatic 
environment  containing  decomposable  organic  matter  is  curtailed,  profound 
changes  occur  both  in  the  composition  of  the  indigenous  microbial  com¬ 
munity  and  in  the  processes  that  these  organisms  carry  out.  The  succes¬ 
sion  of  microorganisms  that  follows  includes  strict  aerobes  first  (that 
require  oxygen)  and  then  facultative  anaerobes  (that  can  exist  with  or 
without  oxygen)  where  both  populations  achieve  their  maximum  levels 
within  the  first  U8  hr  after  flooding  of  a  soil  and  curtailment  of 
oxygen  supply  (Patrick  and  Mikkelsen  1971).  These  are  subsequently 
replaced  by  populations  of  strict  anaerobic  bacteria,  which  cannot  live 
in  the  presence  of  oxygen,  once  the  supply  of  oxygen  within  the  system 
has  been  totally  exhausted  (Takai  and  Kamura  1966;  Takai,  Koyama,  and 
Kamura  1956) . 

20.  The  electrochemical  transformations  that  occur  in  flooded 
soils  and  sediments  have  been  well  documented.  The  oxidation-reduction 
potential  declines  from  highly  positive  preflood  levels  immediately 
after  waterlogging  as  a  consequence  of  the  activities  of  aerobic  and 
facultatively  anaerobic  bacteria.  The  drop  in  oxidation-reduction  poten¬ 
tial  to  negative  values  then  permits  the  initiation  of  full  activities 

by  populations  of  anaerobic  and  sulfate-reducing  bacteria  (Takai  and 
Kamura  1966).  Most  oxidation-reduction  processes  encountered  in  natural 
systems  require  biological  mediation  (Morris  and  Stumm  1967),  which  means 
that  approaches  to  physicochemical  equilibrium  depend  strongly  on  micro¬ 
bial  activities. 

21.  Wetzel  (1975)  has  described  in  considerable  detail  the  many 
changes  in  organic  matter  that  occur  in  association  with  the  microbially 
mediated  processes  of  anaerobic  decomposition  in  sediments.  In  the  sedi¬ 
ments  of  freshwater  systems  such  as  streams,  lakes,  and  ponds,  methano- 
genesis  is  the  major  process  by  which  the  removal  of  large  quantities  of 
carbon  occurs.  The  breakdown  of  organic  matter  to  carbon  dioxide  and 
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methane  occurs  as  a  two-step  process,  although  no  distinction  between 
the  steps  occurs  in  nature  and  the  steps  proceed  simultaneously.  The 
first  step  involves  the  hydrolysis  and  fermentation  of  carbohydrates, 
fats,  and  proteins  to  yield  organic  acids,  alcohols,  and  a  variety  of 
inorganic  products  (Deyl  1961,  McCarty  1964,  Sokatch  1969,  and  Thimann 
1963).  The  organic  acids  and  alcohols  are  then  subject  to  decomposi¬ 
tion  by  the  strictly  anaerobic  methanogenic  bacteria  to  yield  methane 
and/or  carbon  dioxide  (Mah  et  al.  1977)-  The  process  is  completed  when 
gaseous  methane  bubbles  up  out  of  the  sediments  (ebullition)  and  rises 
through  the  water  column,  often  creating  "methane  convection"  currents 
that  may  carry  sediment  particles  and  dissolved  materials  of  the  hypo- 
limnion  with  them  into  the  trophogenic  zone  of  a  lake  or  pond  (Wetzel 
1975)*  Methane  may  also  be  oxidized  by  the  activities  of  methane- 
oxidizing  bacteria  in  portions  of  the  water  column  having  trace  quanti¬ 
ties  of  oxygen  (Wetzel  1975;  Reeburgh  and  Heggie  1977). 

Chemical  effects 

22.  Takai,  Koyama,  and  Kamura  (1956)  incubated  freshly  flooded 
soils  and  found  that  NO^-N  disappeared  at  the  beginning  of  incubation; 
this  was  followed  by  iron  (Fe)  reduction  and  then  active  sulfide  forma¬ 
tion.  Takai  and  Kamura  (1966)  subsequently  found  that  Fe  reduction  did 
not  commence  until  manganese  (Mn)  reduction  was  essentially  complete. 
Turner  and  Patrick  (1968)  found  that  reduction  in  soils  was  not  strictly 
sequential  because  Mn  reduction  was  not  complete  when  Fe  reduction  began. 
However,  little  reduction  of  Mn  and  almost  no  reduction  of  Fe  occurred 
until  both  oxygen  and  nitrate  had  disappeared.  Ponnamperuma  (1972), 
summarizing  the  results  of  many  studies,  reported  that  reduction  of  a 
submerged  soil  proceeds  roughly  in  the  sequence  predicted  by  thermo¬ 
dynamics.  Oxygen  is  first  reduced,  followed  by  nitrate  reduction  after 
oxygen  concentration  has  dropped  to  levels  below  0.2  mg/Jl.  Manganese 
dioxide  is  reduced  next  and  followed  by  reduction  of  iron  and  sulfate. 
Ponnamperuma  (1972)  presented  a  table  that  stated  the  critical  potentials 
for  reduction  of  inorganic  electron  acceptors  in  stirred  soil  suspen¬ 
sions.  These  critical  potentials,  taken  from  the  work  of  Patrick  (1964), 
Connell  and  Patrick  (1968),  and  Turner  and  Patrick  (1968),  are  presented 
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in  the  following  tabulation  and  provide  a  rough  guide  to  the  sequence 
of  reduction  steps  in  natural  systems. 


Observation* 

Eh,7  mV** 

Oxygen  (undetectable) 

330 

Nitrate  (undetectable) 

220 

Reduced  manganese 
(detectable) 

200 

Reduced  iron 
(detectable) 

-120 

Sulfate  (undetectable) 

-150 

Note:  Information  modified  from 

Ponnamperuma  (1972). 

*  Comments  in  parentheses  indicate 

major  changes  in  the  observations 
that  become  evident  at  the  Eh? 
indicated. 

**  Eh,7  mV  =  oxidation-reduction  po¬ 
tential  in  millivolts  at  pH  7* 

23-  Other  results  indicate  that,  even  in  complex  natural  systems 
such  as  flooded  soils,  the  thermodynamic  properties  of  chemical  constit¬ 
uents  apparently  influence  the  sequence  of  reduction.  The  sequence  of 
reduction  in  any  natural  system  should  therefore  be  at  least  qualita¬ 
tively  predictable. 

Measurement  of 

oxidation-reduction  potential 

24.  The  measurement  of  oxidation-reduction  potential  in  natural 
systems  is  replete  with  significant  potential  sources  of  error,  and  good 
technique  is  essential  for  obtaining  reliable  data.  Morris  and  Stumm 
(1967)  pointed  out  that  microenvironments  differing  from  the  overall 
environment  may  be  formed  in  natural  systems.  Stumm  and  Morgan  (1970 ) 
emphasized  that  natural  waters  are  in  a  highly  dynamic  state  rather  than 
in  or  near  equilibrium.  Consequently,  no  single  oxidation-reduction 
electrochemical  potential  (Eh)  representative  of  any  one  of  the  possible 
systems  present  can  be  expected  to  be  reliably  established  with  the 
measuring  electrode  (normally  platinum).  Furthermore,  it  is  generally 
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recognized  that  Eh  measurements  in  oxygenated  natural  waters  are  invalid 
because  of  the  instability  of  the  platinum  electrode  in  the  presence  of 
molecular  oxygen  and  the  dominance  of  oxygen  as  an  electron  acceptor  in 
oxidation-reduction  reactions,  Bohn  (1971)  felt  that  natural  systems 
rarely  achieve  oxidation-reduction  equilibrium  because  of  the  continual 
addition  of  electron  donors  (organic  matter). 

25-  Despite  the  difficulties  involved  in  the  determination  and 
interpretation  of  Eh  readings  in  natural  media,  it  is  widely  accepted 
that  Eh  readings  in  natural  anaerobic  systems  can  provide  valuable  en¬ 
vironmental  information  relative  to  the  dominant  oxidation-reduction 
sequence  in  a  qualitative  sense.  Breck  (197M  considered  t.«ab  measure¬ 
ments  of  Eh  with  a  platinum  electrode  might  depict  trends  in  a  useful 
and  quite  reliable  manner.  Ponnamperuma  (1972)  considered  Eh  a  useful 
environmental  parameter  in  spite  of  its  limitations  ana  cited  numerous 
studies  in  support  of  his  contention.  Whitfield  (  '909;  ■~nd  Eh  useful 
as  a  semiquantitative  indicator  of  the  degree  of  stagnation  of  a  particu¬ 
lar  aquatic  environment  as  did  Ponnamperuma  (1972)  for  flooded  soils  and 
sediments. 


Chemical  oxidation- 
reduction  sequences  in 
natural  aquatic  environments 

26.  The  following  section  examines  the  various  oxidation-reduction 
chemical  systems  present  in  natural  aquatic  environments  in  sequence 
from  most  oxidized  to  most  reduced. 

27.  Oxygen.  The  presence  of  0^  in  a  system  will  maintain  Eh  at 
potentials  in  excess  of  +300  to  +k00  mV.  Platinum  electrodes  do  not, 
however,  register  steady  or  reproducible  redox  potentials  in  poorly 
poised  (unstable)  systems,  such  as  aerated  soils  and  natural  waters 
(Ponnamperuma  1972).  This  problem  in  Eh  measurements  in  aerobic  systems 
has  been  attributed  to  the  absence  of  electroactive  systems  in  suffi¬ 
ciently  high  concentrations  and  to  the  presence  of  mixed  systems  that 
are  not  in  equilibrium  among  themselves  (Morris  and  Stumm  1967;  Stumm 
and  Morgan  1970).  Patrick  and  Mikkelsen  (1971)  reported  that  aerated 
soils  have  characteristic  redox  potentials  in  the  range  of  +U00  to 


+700  mV.  Dissolved  oxygen  does  not  exert  the  potential  that  would  be 
expected  if  it  were  functioning  at  equilibrium  (Garrels  and  Christ  1956). 

Instead,  it  acts  like  a  much  weaker  oxidizing  agent,  with  an  Eh  of  +650 
to  +700  mV  at  pH  =  0  and  about  +300  to  +350  mV  at  pH  =  8  .  Systems 
exposed  to  air  show  potentials  (in  volts  at  25 °C)  according  to  the 

approximate  relation  (Garrels  and  Christ  1956)  Eh  =  0.70  -  0.06  pH  ,  X 

in  contrast  to  the  expected  theoretical  relation  Eh  =  1.23  -  0.06  pH  . 

28.  Thus,  it  is  apparent  that  Eh  measurements  in  aerated  systems 
are  of  little  value.  The  oxidation  state  of  aerobic  systems  is  better 
characterized  by  measuring  the  DO  concentration  rather  than  by  determin¬ 
ing  the  Eh. 

29-  Nitrogen.  It  has  been  established  that  the  Eh  of  anaerobic 
systems  will  be  stabilized  at  potentials  between  +200  and  +400  mV  as 
long  as  NO^  is  present  (Turner  and  Patrick  1968,  Ponnamperuma  1965, 

Ponnamperuma  and  Castro  1964,  Yamane  and  Sato  1968,  and  Yamane  1969)* 

_  +2 
The  presence  of  N0_  prevents  the  production  of  other  reduced  forms  (MN  , 

+2  —2 

FE  ,  S  ,  CH^,  and  E^)  in  the  sequence  of  oxidation-reduction  events  in 
lake  waters  and  muds  (Pearsall  and  Mortimer  1939;  Mortimer  1941  and  1942) 
and  in  submerged  soils  (Turner  and  Patrick  1968;  Ponnamperuma  1972). 

Consequently,  the  presence  of  N0~  will  poise  the  system  at  a  slightly 
reduced,  but  relatively  oxygen-free  level. 

30.  The  NO^-NOg  couple  has  a  calculated  Eq  of  430  mV  at  pH  7 
(Whitfield  1972).  Many  workers  (Bohn  1971,  Stumm  1967,  and  Whitfield 
1972),  however,  consider  the  NO^-NO^-NH^  couples  irreversible  since  the 
reduction  of  NQ^  to  NO^  to  NH^  and/or  N^O  and  Ng  (denitrification)  is 
normally  mediated  by  microorganisms.  Evidence  exists  that  although  N0~ 
is  subject  to  irreversible  reactions  such  as  denitrification,  it  will 
increase  the  Eh  and  retard  a  fall  in  Eh  (Engler  and  Patrick  1974)  when 
NO^  is  added  to  anaerobic  systems.  Nitrate  may  affect  the  Eh  by  means 
of  an  electroactive  mediator  system.  Elema  (1932),  cited  by  Thimann 
(1963),*  found  that  a  steady  Eh  was  established  when  solutions  of  nitrate 

*  Elema,  B.  1932.  De  bepaling  von  oxydatie-reductie  potentiaol  in 
bacteriencultures  en  here  beteekenis  voor  de  stof ivisseling,  Ph.  D. 

Dissertation,  Delft,  The  Netherlands. 
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or  a  nitrite-containing  hydrogen  donor  (usually  ethanol)  were  inoculated 
with  a  culture  of  denitrifying  microorganisms.  Elema  concluded  that  the 
reversible  substance  causing  a  steady  Eh  must  have  been  released  from 
the  bacteria  because  the  nitrite  and  alcohol  alone  were  unable  to  regis¬ 
ter  a  steady  potential.  Moreover,  enzymes  (such  as  those  in  bacteria) 
are  known  to  have  the  capacity  to  mediate  the  potentials  of  redox  couples 
that  do  not  exhibit  activity  by  themselves  (Clark  1972). 

31.  Manganese.  Manganese  is  not  as  effective  as  nitrate  in  stabi- 

+4 

lizing  or  poising  Eh.  Once  Mn  reduction  had  commenced,  the  influence 
of  MnO^  in  stabilizing  the  Eh  was  found  to  be  weak  because  of  its  rela¬ 
tive  insolubility  and  because  of  its  use  as  an  electron  acceptor  in 
respiration  by  only  a  limited  number  of  bacteria  (Ponnamperum  1972). 

Nhung  and  Ponnamperuma  (1966)  found,  however,  that  added  MnO^  was  much 
more  effective  than  added  Fe(OH)^  in  retarding  the  fall  in  Eh  of  a 
flooded  acid  sulfate  soil.  The  effectiveness  of  the  MnO^  in  retarding 
soil  reduction  was  credited  to  the  high  standard  oxidation-reduction 
potential  (Eq  =  1.230  mV)  of  the  MnO^  Mn+2  system.  Native  or  added  MnO,,, 
in  addition  to  retarding  the  decrease  in  Eh  or  the  onset  of  moderately 

anaerobic  conditions  in  flooded  soils,  was  also  found  to  prevent  the 

+2 

buildup  of  high  concentrations  of  Fe  and  other  reduction  products 

(Ponnamperuma  and  Castro  1964;  Ponnamperuma,  Yaun,  and  Nhung  1965). 

+4 

32.  Many  studies  have  focused  on  the  effect  of  the  Mn  reduction 

process  on  Eh.  Ponnamperuma,  Loy,  and  Tianco  (1969 )  found  that  the  E  's 

o 

of  Mn  compounds  in  soils  and  sand  plus  MnO^  were  markedly  lower  than 
those  of  the  ideal  oxides  and  varied  with  the  soil  and  the  intensity  of 
reduction.  Experimentally  derived  Eq's  deviated  least  from  theory 
and  were  closest  to  those  of  the  sand-MnO^  medium  for  soils  high  in  Mn 
and  low  in  organic  matter.  Ponnamperuma,  Loy,  and  Tianco  (1969)  used 

the  constancy  of  E  in  a  soil  despite  variations  in  Eh,  pH,  PCO  ,  and 

+2 

Mn  activity,  as  the  criterion  for  the  presence  of  a  particular  Mn 

system.  Olomu,  Racy,  and  Cho  (1973)  found  that  E  values  calculated 

+2  ° 

from  experimentally  measured  Eh,  pH,  and  Mn  concentrations  appeared 

+2 

to  indicate  the  presence  of  the  M^O^-Mn  redox  system. 

33.  Other  workers  have  not  found  any  relationship  between  Eh  and 
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the  Mn+VMn+2  couple,  Bohn  (1968)  found  that  the  electrode  potential 
of  two  soils  was  independent  of  the  Mn+1+/Mn+2  couple.  Gotoh  and  Patrick 
(1972)  found  marked  deviation  from  the  equations  calculated  for  pure 

systems  for  the  calculation  of  equations  that  relate  experimentally 

+2 

obtained  Mn  activity  to  pH  and  Eh  in  the  pH  range  of  6  to  8.  Bohn 

+2 

(1970)  found  that  Mn  concentrations  in  control  and  autoclaved  soil 

+2 

suspensions  differed  widely  from  the  Mn  concentrations  calculated  from 
models;  the  latter  require  measurement  of  Eh  or  assume  Eh  to  be  fixed 
by  another  redox  couple.  However,  the  soil  suspensions  studied  (Bohn 
1970)  were  kept  in  loosely  stoppered  containers  and  were  therefore 
subject  to  air  oxidation. 

31*.  Iron.  Takai,  Koyama,  and  Kamura  (1963)  found  a  close  connec¬ 
tion  between  Eh  and  pH  of  flooded  soils  and  the  amount  of  ferrous  iron 
in  the  soil.  Takai  and  Kamura  (1966)  considered  oxidized  iron  in  soils 
to  be  the  most  important  factor  retarding  the  onset  of  highly  reducing 
conditions  because  of  the  high  concentrations  of  reducible  iron  that 

normally  occur  in  most  soils.  Ponnamperuma  (1972)  reported  that  the  Eh, 

+2 

pH,  and  Fe  concentrations  in  reduced  lake  waters  reported  by  various 

workers  (Mortimer  19H1  and  19^2)  conformed  to  the  equilibrium  values  of 

+2 

the  Fe(OH)^  Fe  system.  Hem  (i960)  found  that  the  behavior  of  iron 

in  natural  systems  was  predicted  to  a  large  extent  by  fluctuations  in 

the  oxidation-reduction  potential  and  pH  environment.  Ponnamperuma, 

Tianco,  and  Loy  (1967)  confirmed  the  participation  of  Fe(OH)^  and 

Fe^(OH)^  in  redox  equilibria  in  flooded  soils. 

35.  Ponnamperuma,  Tianco,  and  Loy  (1967)  found  that  the  Eh,  pH, 

+2  .  , 

and  Fe  activity  of  32  flooded  soil  solutions  sampled  weekly  for  17 
weeks  conformed  closely  to  the  equation 

Eh  =  1.058  -  0.059  log  Fe+2  -  0.177  pH  (l) 

during  the  period  of  submergence,  up  to  the  point  where  water-soluble 

+2 

Fe  reached  a  peak  when 

Eh  =  0.L29  -  0.059  log  Fe+2  -  0.236  pH  (2) 
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and  subsequently 


Eh  =  0.U29  -  0.059  pH  (3) 

better  expressed  the  situation.  This  confirmed  the  participation  of 
Fe^H)^  and  Fe^(OH)g  in  redox  equilibria  in  flooded  soils.  The  partici¬ 
pation  of  iron  in  redox  equilibria  of  flooded  soils  was  also  supported 

by  the  work  of  Gotoh  and  Patrick  (197^),  whose  results  suggested  that 

+2 

Fe  activity  in  a  flooded  soil  might  be  governed  generally  by  the  solu¬ 
bility  of  ferric  oxyhydroxide .  Van  Breeman  (1969)  claimed  that  the 
equilibrium  between  ferric  and  ferrous  iron  was  largely  governed  by  ill- 
defined  ferric  oxides,  which  were  intermediate  in  stability  between  amor¬ 
phous  Fe(OH)^  and  a  -  FeOOH  . 

36.  Other  workers  have  been  unable  to  show  the  participation  of 

iron  in  redox  equilibria.  Olomu,  Racy,  and  Cho  (1973)  found  that  the 

+2 

calculated  values  of  Eq  ,  assuming  an  operative  Fe^O^-Fe  redox  system, 
were  far  greater  than  the  theoretical  value.  However,  they  also  found 
that  almost  all  the  Fe  in  soil  solution  was  complexed  with  organic  mat¬ 
ter.  The  high  degree  of  Fe  complexing  in  soil  solution  was  suggested  as 

the  reason  for  the  failure  to  obtain  the  theoretical  E  value  for  the 

o 

Fe^O^-Fe  redox  system  from  the  experimental  data.  Bohn  (1969)  found 

that  iron(ll)  sulfate  added  to  aerobic  soil  suspensions  did  not  govern 

+3 

Eh  because  of  the  immediate  precipitation  of  the  Fe  formed.  Coursier 
(1952)  found  that  the  concentrations  of  iron(Il)  and  (III)  each  must 
exceed  10  ^  M  (O.56  mg/8.)  in  order  to  govern  Eh.  In  addition,  Collins 
and  Buol  (1970)  showed  that  some  Mn  and  Fe  precipitation  reactions  were 
not  reversible  over  short  periods  of  time. 

37-  The  effect  of  the  iron  system  is  not  as  great  as  those  of 

nitrate  and  magnanese  on  Eh.  Ponnamperuma  (1972)  stated  that  because 

+2 

of  the  low  standard  potential  of  the  Fe(0H)^-Fe  system,  its  influence 
was  not  as  obvious  as  that  of  N0~  or  MnO^.  Nhung  and  Ponnamperuma 
(1966)  found  that  added  Fe(OH)^  retarded  to  a  limited  extent  the  fall  of 
Eh  in  a  flooded  acid  sulfate  soil.  The  effect  of  the  added  Fe(OH)^  on 
Eh  was  not  nearly  as  great  as  that  of  MnO^.  Ponnamperuma  (1965)  noted 
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that  soils  high  in  Fe  showed  a  slower  rate  of  fall  in  Eh  in  the  zone 

„  +3 

-50  to  -200  mV  than  did  soils  low  in  Fe 

-2  -2 

38.  Sulfur.  The  reduction  of  SO^  to  S  does  not  occur  until  an 
Eh  of  approximately  -150  mV  is  reached  within  a  pH  range  of  6.8  and  7.1 
(Connell  and  Patrick  1968)  and  represents  intense  reducing  conditions. 
Martens  and  Berner  (197M  have  recently  shown  that  methane  concentrations 
in  the  interstitial  waters  of  Long  Island  Sound  sediments  did  not  reach 
appreciable  concentrations  until  about  90  percent  of  seawater  sulfate 
was  removed  by  sulfate-reducing  bacteria.  Laboratory  studies  of  anoxic 
marine  sediments  indicated  that  methane  production  does  not  occur  until 
dissolved  sulfate  is  totally  exhausted  (Martens  and  Berner  197M* 

39.  Berner  (1963)  demonstrated  that  the  Eh  of  many  sediments 
containing  H^S  is  controlled  by  the  reversible  half-cell 

HS"  S°  ,  +  H+  +  2e~  (1) 

aq  rhmb  aq 

which  is  equivalent  to  the  half-cell 


s  v,  v  +  2e 
rhmb 


The  controlling  S“  /S°,  ,  half-cell  found  by  Berner  (1963)  was  truly 
aq  rhmb 

reversible  (Voge  19^9;  Peschanski  and  Valenski  19^9).  Other  workers 
(Whitfield  1969;  Skopintsev,  Romenskaya,  and  Smirnov  1966)  have  also 
demonstrated  that  the  platinum  electrode  responds  to  the  sulfide-sulfur 
half-cell  in  a  wide  variety  of  reduced  systems  containing  sulfur. 

40.  The  Eh  of  a  natural  H^S  containing  sediment  is  controlled  by 
the  sulfide-sulfur  half-cell  only  if  an  amount  of  oxidizing  agent  suffi¬ 
cient  to  bring  about  the  oxidation  of  sulfide  to  elemental  sulfur  is 
present  (Berner  1963).  The  presence  of  ferric  hydroxide  or  geothite  will 
oxidize  HgS  to  elemental  sulfur;  from  there  further  oxidation  to  a 
saturated  polysulfide  solution  can  occur  (Berner  1963). 

Ul.  Hydrogen.  Reducing  conditions  become  sufficiently  strong  in 
some  anaerobic  environments  for  the  H+-H2  couple  to  affect  the  Eh. 

Yamane  and  Sato  (1968)  have  found  that  the  first  minimum  potential 


following  flooding  can  be  as  low  as  -U20  mV  and  can  be  accompanied  by 
the  evolution  of  hydrogen.  Baas-Becking,  Kaplan,  and  Moore  (i960) 
reported  that  the  H+-H2  exchange  current  was  great  enough  for  its  poten¬ 
tial  to  be  the  major  determinant  of  Eh  in  strongly  reducing  environments 
The  influence  of  the  couple  on  Eh  at  25°C  may  be  described  by  the 

equation 


Eh  =  -0.029  pH2  -  0.059  pH  (6) 

U2.  The  hydrogen  ion  concentration,  in  addition  to  the  H+-H2 
couple,  may  also  influence  Eh.  Clark  (1972)  reported  that  the  hydrogen 
ion  concentration  affects  Eh  by  influencing  the  dissociation  of  oxidant 
or  reductant.  The  effect  of  pH  on  Eh  has  been  reviewed  in  detail  by 
Laitinen  (i960)  and  Ponnamperuma  (1972). 

Effect  of  organic  mat- 
ter  on  anaerobic  processes 

U3.  Organic  matter,  in  addition  to  serving  as  an  energy  source 
for  bacteria,  is  also  a  source  of  electrons  (electron  donor)  for  biologi 
cal  reduction  processes  (Ponnamperuma  1972).  The  organic  matter  that 
serves  as  a  source  of  electrons  becomes  oxidized  and  can  potentially 
form  a  couple  with  its  reduced  precursor.  Thus,  it  is  not  surprising 
that  the  addition  of  decomposable  organic  matter  will  affect  the  rate 
and  intensity  of  reduction  in  anaerobic  systems.  Addition  of  organic 
matter  to  anaerobic  systems  will  generally  increase  the  rate  of  reduc¬ 
tion  (Ponnamperuma  1972;  Takai,  Koyama,  and  Kamura  1963).  A  rapid  de¬ 
cline  in  Eh  was  characteristic  of  soils  with  low  contents  of  reducible 
Fe  and  Mn  and  a  high  organic  matter  content  (Patrick  and  Mahapatra  1968) 
Ponnamperuma  and  Castro  (I961*)  suggested  that  the  potentials  of  reduced 
soils  might  be  potentials  of  organic  systems  mediated  by  enzymes  and 
bacteria. 

U4.  Organic  matter  can  exert  a  decided  impact  on  the  Eh  of  a 
system  despite  the  electrochemical  irreversibility  of  most  organic 
couples  (Clark  1972),  since  irreversible  redox  couples  can  exert  an 
electrode  potential  via  a  mediator.  Lehmann  (1930),  reviewed  by  Clark 


(1972),  demonstrated  that  the  irreversible  suceinate/fumarate  system 
could  influence  electrode  potentials  in  the  presence  of  a  mediator. 
Electrode  potentials  at  constant  pH  conformed  to  Eh's  predicted  from 
the  ratio  of  [ succinate ]/[fumarate]  and  were  independent  of  the  nature 
or  concentration  of  the  mediator.  Within  a  limited  range  of  pH,  the 
potentials  were  in  conformity  with  the  theoretical  slope  of  0.06. 

45 .  When  a  soil  is  submerged,  soil  solution  concentrations  of 
volatile  organic  acids  increase,  reach  a  peak  of  value  of  10  to  1+0  mmol/Jl 
in  1  to  2  weeks,  and  then  decline  to  less  than  1  mmol/Jl  a  few  weeks 
later  ( Ponnamperuma  1972).  Low  temperatures  lead  to  extensive  accumula¬ 
tion  of  organic  acids  in  the  soil  solution  (International  Rice  Research 
Institute  (IRRI)  1969).  Motomura  (1962),  as  cited  by  Ponnamperuma  (1972), 
found  that  soils  high  in  native  or  added  organic  matter  have  correspond¬ 
ingly  high  concentrations  of  organic  acids.  There  is  evidence  that  this 
buildup  of  organic  acids  may  affect  Eh.  Clark  (1972)  summarized  results 
of  many  studies  that  showed  that  addition  of  the  alkali  salts  of  acetic, 
citric,  lactic,  malonic,  oxalic,  salicyclic,  or  tartaric  acid  generally 
decreased  the  Eh  of  an  Fe  /FE  system  at  constant  pH.  These  observa¬ 
tions,  coupled  with  the  buildup  of  organic  acids  in  anaerobic  systems, 
further  demonstrate  the  potential  effects  that  organics  may  exert  on 

Eh. 

Discussion 

46.  The  survey  of  the  literature,  coupled  with  other  studies 
(Brannon  et  al.  1978),  suggests  that  the  microorganisms  and  anaerobic 
processes  involved  in  the  decomposition  of  organic  matter  and  in  the 
reduction  of  nitrate,  manganese,  iron,  and  sulfate  in  organic-laden 
waters  are  very  similar  to  those  occurring  in  flooded  soils  and  sedi¬ 
ments.  The  sequence  of  reduction  of  these  inorganic  electron  acceptors 
generally  corresponded  to  results  obtained  by  other  workers  (Morris  and 
Stumm  1967;  Takai  and  Kamura  1966;  Takai,  Koyama  and  Kamura  1956; 

Connell  and  Patrick  1968;  Patrick  1964;  Ponnamperuma  1972;  Turner  and 
Patrick  1968)  wherein  dissolved  oxygen  was  rapidly  depleted  in  the 
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culture  medium  with  the  rate  of  removal  being  proportional  to  the  amount 
of  organic  matter  present  in  the  system. 

U7.  Previous  studies  by  the  authors  of  this  report  (Brannon  et  al. 
1978)  have  shown  that  the  succession  of  microorganisms  occurring  in 
aquatic  microcosms  closely  follows  the  successional  pattern  observed  by 
other  workers  in  flooded  soils  and  sediments.  The  succession  of  micro¬ 
organisms  in  flooded  soils  includes  aerobes  first  and  facultative 
anaerobes  second,  where  both  populations  achieve  their  maximum  levels 
within  the  first  48  hr  after  flooding  (Ponnamperuma  1972).  During  this 
period,  the  supply  of  oxygen  in  the  sediment  is  totally  exhausted,  and 
populations  increase  of  strict  anaerobic  bacteria,  to  which  0^  is  a 
lethal  toxin  (Takai  and  Kamura  1966;  Takai,  Koyama,  and  Kamura  1956). 

48.  Nitrate  reduction  commences  when  0^  is  either  partially  or 
totally  consumed,  with  denitrification  probably  serving  as  the  main 
removal  mechanism  for  nitrate. 

49.  Brannon  et  al.  (1978)  has  clearly  demonstrated  the  similarity 
between  reduction  processes  occurring  in  systems  simulating  reservoir 
waters  and  those  occurring  in  soils  and  sediments,  indicating  that  re¬ 
sults  of  studies  in  flooded  soils  may  well  be  applicable  to  other  sys¬ 
tems  involving  lake  sediments.  Indeed,  Ponnamperuma  (1972)  has  studied 
the  movement  of  the  boundary  layer  between  the  aerobic  and  anaerobic 
zones  from  within  lake  sediments  to  the  water  above  the  mud  surface ,  and 
into  the  hypolimnion.  By  this  process,  the  water  of  the  hypolimnion  may 
become  an  extension  of  the  sediment  interstitial  water,  although  the 
intensity  of  reduction  and  the  reduction  products  observed  may  not  be 
identical.  Mortimer  (1941,  1942)  has  observed  that,  in  thermally 
stratified  eutrophic  lakes,  the  Eh  declines  rapidly  in  the  thermocline, 
and  in  the  hypolimnion  the  Eh  may  acquire  values  as  low  as  those  in  the 
underlying  reduced  sediments. 

50.  Results  of  the  Brannon  et  al.  (1978)  studies  indicate  that 

the  amount  of  organic  matter  present  will  be  a  major  determinant  of  the 

+3  +4 

rate  of  reduction  of  Fe  and  Mn  .  The  concentrations  that  reduced 
components  reach  in  the  system  depend  not  only  on  the  amount  of  organic 
matter  present,  but  also  on  the  amount  of  reducible  substrate  initially 
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present.  For  accurate  predictions  of  anaerobic  processes  in  hypolimnions 

of  CE  reservoirs,  the  amount  of  decomposable  organic  matter  and  reducible 
+3  +L 

Fe  and  Mn  must  be  determined. 

51.  To  illustrate  the  sequential  nature  of  reductive  processes 
as  they  occur  in  natural  systems,  a  series  of  seven  figures  (Figures  3-9) 
is  presented  to  illustrate  aerobic/anaerobic  processes  observed  in 
waters  overlying  the  soil  horizon  taken  from  the  area  of  the  proposed 
Twin  Valley  Lake  in  western  Minnesota  (all  figures  are  taken  from 
Gunnison  et  al.  (1979)).  In  each  of  the  figures,  depletion  of  DO  and 
establishment  of  anaerobic  conditions  occur  in  the  same  order  as  indi¬ 
cated  in  Figure  3;  the  "Transition"  section  marked  between  lU  and  19  days 
of  incubation  in  each  of  the  figures  denotes  a  changeover  from  aerobic 
to  anaerobic  conditions.  Thus,  as  oxygen  is  depleted  from  the  system 
and  anaerobic  conditions  ensue,  pH  declines  from  9.*+  to  6.8  while  conduc¬ 
tivity  increases  (Figure  L);  total  dissolved  organic  carbon  (TOC)  accumu¬ 
lates  in  the  water  in  an  erratic  manner  while  total  inorganic  carbon  (TIC) 


Figure  3-  Changes  in  DO  concentration  and  oxidation-reduction  (redox) 
potential  in  water  columns  of  the  reactor  units  during  100-day  incuba¬ 
tion  period.  Bars  around  each  mean  value  represent  standard  error  of 

the  mean 


Figure  h.  Changes  in  pH  and  conductivity  in  the  water  columns  of  the 
reactor  units  during  100-day  incubation  period.  Deviations  from  mean 
values  were  insignificant  and,  therefore,  are  not  presented  here 

increases  nearly  linearly  (Figure  5);  nitrate-nitrogen  (NO^-N)  is  de¬ 
pleted  steadily  as  anaerobic  conditions  are  developing;  nitrite-nitrogen 
(NO^-N)  appears  in  a  transitory  fashion  prior  to  the  onset  of  anaerobic 
conditions,  and  ammonium  nitrogen  (NH^-N)  does  not  accumulate  until  after 
anaerobic  conditions  have  been  established  (Figure  6);  total  phosphorus 
accumulates  in  the  water  column  in  an  irregular  manner,  but  orthophos¬ 
phate  phosphorus  does  not  accumulate  in  significant  concentrations  until 
after  the  establishment  of  anaerobic  conditions  (Figure  7)  and  iron 
and  manganese  release  from  the  soil  to  the  water.  Each  of  the  latter 
constituents  requires  anaerobic  conditions  for  accumulation  (Figure  8), 
while  sulfide  production  from  sulfate  does  not  ensue  until  anaerobic 
conditions  have  become  well  established  (Figure  9)* 


CONCENTRATION  OF 
DISSOLVED  ORGANIC  CARBON,  m  ft 


Figure  7.  Changes  in  total  P  and  ortho-P  in  the  water  columns  of  the 
reactor  units  during  100-day  incubation  period.  Bars  around  each 
mean  value  represent  standard  error  of  the  mean 


Figure  8.  Changes  in  concentration  of  soluble  reduced  manganese 
(total  Mn)  and  iron  (total  Fe)  in  the  water  columns  of  the  reac¬ 
tor  units  during  100-day  incubation  period.  Bars  around  each 
mean  value  represent  standard  error  of  the  mean 
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INCUBATION  PERIOD,  DAYS 

Figure  9-  Changes  in  concentration  of  dissolved  sulfate  and  total 
suspended  sulfide  in  the  water  columns  of  the  reactor  units  during 
100-day  incubation  period.  Bars  around  each  mean  value  represent 

standard  error  of  the  mean 
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PART  IV:  DOSAP-DISSOLVED  OXYGEN  SUBROUTINE 
WITH  ANAEROBIC  PROCESSES 


52.  The  DOSAP  formulation  described  herein  was  developed  as  the 
conceptual  basis  for  a  mathematical  subroutine  for  simulating  anaerobic 
processes.  Current  mathematical  water  quality  models  lack  the  capability 
to  track  anaerobic  processes.  Although  intended  for  use  with  the  CE 
reservoir  water  quality  model,  CE-QUAL-R1,  DOSAP  can  be  interfaced  with 
other  numerical  models. 

53.  This  portion  of  the  report  describes  in  two  sections  the 
basic  structure  of  the  DOSAP  formulation.  The  first  section  describes 
the  overall  organization  of  the  subroutine,  its  various  components,  and 
then  integrates  these  components  into  a  functionally  unified  picture. 

The  second  section  analyzes  each  component  in  detail  and  explains  the 
various  physical-chemical  subsections  of  each  component.  A  compilation 
of  representative  rate  data  for  the  subroutine,  as  obtained  from  the 
literature  and  the  authors'  studies,  is  presented  and  discussed  in  Part  V. 

Subroutine  Organization  and  Function 

51*.  The  DOSAP  subroutine  consists  of  a  series  of  seven  steps  or 
phases  that  gradually  and  continuously  proceeds  from  aerobic  (oxic)  to 
anaerobic  (anoxic)  conditions  and  then  proceeds  through  a  limited  suc¬ 
cession  of  anerobic  product-generating  loops ,  each  of  which  is  more  in¬ 
tensely  anaerobic  (reducing)  than  its  predecessors. 

Phases  of  DOSAP 

55-  Phase  I:  DO  depletion.  This  phase  is  initiated  at  the  onset 
of  thermal  stratification  in  the  reservoir.  During  stratification,  the 
diffusive  exchange  of  DO  between  the  epilimnion  and  the  hypolimnion  and 
advective  transport  of  oxygen-bearing  water  into  the  hypolimnion  are 
greatly  diminished.  Decomposition  of  organic  matter  proceeds  together 
with  the  oxidation  of  reduced  anaerobic  constituents  released  from  reser¬ 
voir  bottom  sediments.  Through  these  processes,  water  column  DO  may  be 
depleted  down  to  a  level  of  2  mg/ l,  the  initiation  point  for  Phase  II. 


56.  Phase  II:  DO  exhaustion.  Depletion  of  DO  continues  in  this 
phase  until  all  of  the  DO  has  been  consumed;  this  is  the  signal  to 
initiate  Phase  III.  As  a  consequence  of  the  decreasing  concentration  of 
DO,  nitrate  begins  to  replace  oxygen  as  an  inorganic  electron  acceptor 
for  microbial  respiration  processes.  As  a  result,  the  concentration  of 
nitrate  begins  to  decrease  during  Phase  II. 

57*  Phase  III:  Ammonium  accumulation.  Nitrate  is  rapidly  de¬ 
pleted  and  lost  to  the  system.  Ammonium  begins  to  accumulate  in  the 
water  column.  This  is  a  consequence  of  the  nonavailability  of  oxygen 
to  be  used  in  the  microbial  oxidation  of  ammonium  that  is  derived  from 
decomposition  of  nitrogenous  organic  matter  and  from  anaerobic  reduction 
of  nitrate  to  ammonium.  Once  all  nitrate  has  been  depleted,  the  system 
shifts  to  Phase  IV. 

58.  Phase  IV:  Manganese  accumulation.  Soluble  reduced  manganese 
released  from  the  sediments  begins  to  accumulate  in  the  hypolimnetic 
water  column.  The  oxidation-reduction  sequential  process  continues 
until  reducing  conditions  become  more  intense,  as  indicated  by  maximum 
rates  of  manganese  reduction  and  accumulation.  A  decrease  in  the  rate 
of  manganese  reduction  signals  the  start  of  Phase  V. 

59-  Phase  V:  Iron  accumulation.  Reduction  of  ferric  iron  (Fe  ) 

2+ 

to  ferrous  iron  (Fe  )  begins  and  the  subsequent  accumulation  of  soluble 
ferrous  iron  in  the  water  column  continues  until  the  rate  of  iron  accumu¬ 
lation  begins  to  approach  zero.  Phase  VI  is  then  initiated. 

60.  Phase  VI:  Sulfate  reduction.  Reduction  of  sulfate  present 
in  the  water  and/or  any  sulfate  entering  the  hypolimnion  with  inflows  is 
initiated,  and  sulfide  begins  to  appear  in  the  water  column.  Sulfide 
released  from  the  sediment  to  the  hypolimnion  is  stable  due  to  the  in¬ 
tense  anaerobic  conditions.  The  rotten  egg  odor  characteristic  of  hy¬ 
drogen  sulfide  will  be  evident  in  any  hypolimnetic  waters  that  reach  the 
surface,  either  upon  destratification  and  mixing  of  the  reservoir,  or 
in  samples  of  the  hypolimnion  that  are  brought  up  for  examination.  In 
addition,  sulfide  and  its  attendant  odors  will  also  occur  in  any  releases 
made  from  the  hypolimnion.  Finally,  the  sulfide  will  easily  combine  with 
any  reduced  iron  present  in  the  water  to  form  insoluble  iron  sulfide, 
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which  is  then  precipitated.  This  phase  is  terminated  when  the  supply  of 
sulfate  has  been  exhausted,  and  Phase  VII  is  initiated. 

61.  Phase  VII:  Methanogenesis.  Methane  formation  is  initiated 
and  continues  until  the  system  is  perturbed  or  the  reservoir  destratifies . 

62.  Destratification.  As  an  interim  measure,  the  destratification 
process  automatically  reoxidizes  all  reduced  chemical  entities  dissolved 
or  suspended  in  the  layer.  While  this  simplified  approach  does  not  pre¬ 
vent  the  useful  application  of  the  subroutine,  certain  inconsistencies, 
such  as  the  instantaneous  oxidation  of  manganese,  result  because  of  the 
automatic  reoxidation  process.  However,  the  presence  of  such  inconsis¬ 
tencies  is  minimal  relative  to  the  error  introduced  by  omitting  the 
reoxidation  step.  A  more  realistic  formulation  for  processes  whereby 
reduced  substances  accumulated  under  anaerobic  conditions  are  chemically 
and  biologically  oxidized  after  destratification  is  presently  being 
developed  and  tested  as  part  of  this  work  unit.  Once  this  has  been 
completed,  various  interactions  between  dissolved  and  particulate  matter 
will  also  be  added  to  the  subroutine. 


DOSAP  function  in 
relation  to  CE-QUAL-R1 

65.  An  example  of  the  manner  in  which  DOSAP  might  function  as  an 
algorithm  in  CE-QUAL-R1  or  an  equivalent  water  quality  model  is  depicted 
in  Figure  10.  The  hypolimnion  of  the  reservoir  may  be  modeled  as  a 
series  of  horizontal  slices  of  variable  thickness  (user  selected)  stacked 
from  the  bottom  of  the  water  column.  The  DOSAP  algorithm  would  be 
initiated  following  the  onset  of  thermal  stratification.  According  to 
this  example,  all  slices  of  the  hypolimnion  would  be  in  Phase  I  of  the 
DOSAP  algorithm  at  this  time  (Figure  10a). 

6U.  After  a  period  of  thermal  stratification,  the  oxygen  demand 
of  the  bottom  sediments  and  bottom  layers  of  water  may  become  large 
enough  to  deplete  these  layers  of  DO,  thereby  pushing  these  sediments 
and  layers  into  Phase  II.  With  continued  oxygen  depletion,  resulting 
from  oxygen  depletion  exceeding  oxygen  gains,  the  upper  layers  of  the 
hypolimnion  could  also  become  depleted  and  these  layers  would  enter 
Phase  II  as  shown  in  Figure  10b. 


c.  Late  after  onset  of  stratification.  DO  is  gone  and  ammonium, 
manganese,  and  iron  are  being  accumulated  in  the  hypolimnion 

Figure  10.  Diagram  depicting  hypothetical  operations  of  DOSAP  in 
simulating  conditions  in  the  hypolimnion  of  a  reservoir  being 
modeled  as  a  series  of  horizontal  layers.  Stippled  areas  in 
layers  E  through  L  represent  suspended  sediments.  For  details, 

see  text 
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65.  Continued  thermal  stratification  and  development  of  anaerobic 
conditions  would  gradually  result  in  each  of  the  layers  in  the  hypolim- 
nion  progressing  through  the  seven  DOSAP  phases  with  eventual  develop¬ 
ment  of  an  anaerobic  hypolimnion  (Figure  10c).  Hypolimnetic  discharges 
would  contain  dissolved  manganese,  ammonium,  and  traces  of  iron.  Be¬ 
cause  of  the  prolonged  period  of  oxygen  depletion  and  the  large  require¬ 
ment  for  oxidizable  carbon  sources,  it  is  doubtful  that,  in  most  reser¬ 
voirs  ,  the  topmost  layer  of  the  hypolimnion  would  procede  much  beyond 
Phase  III  or  IV.  In  addition,  it  is  doubtful  also  that  Phase  VII  would 
ever  be  found  above  the  bottommost  layers. 

66.  The  DOSAP  algorithm  can  also  be  used  to  simulate  the  effects 
of  a  metalimnetic  oxygen  demand  where  high  oxygen  demand  associated  with 
interflows  of  decomposition  of  settling  organic  particulates  causes 
depletion  of  oxygen  in  a  layer  of  the  metalimnion;  therefore,  this 
layer(s)  enters  Phase  II. 


Phase  Components  and  Interrelationships 


67.  Detailed  descriptions  of  each  DOSAP  phase  are  presented 
below.  The  description  for  each  phase  contains:  (l)  the  pH  range  over 
which  the  phase  is  functional;  (2)  the  oxidation-reduction  (redox)  po¬ 
tential  range  for  the  phase;  (3)  the  DO  concentration,  if  any;  ( U )  a 
synopsis  or  scenario  that  summarizes  the  sequence  of  major  events  that 
occur  during  the  phase;  and  (5)  a  diagram  that  depicts  the  major  chem¬ 
ical  components  and  interactions  present  during  the  phase  along  with  a 
tabulation  of  the  processes  represented  by  each  of  the  arrows  in  the 
diagram  and  a  list  of  algorithms  summarizing  the  various  processes  oper¬ 


ative  on  each  component  in  the  diagram.  The  boxes  surrounding  each  com¬ 
ponent  indicate  whether  that  component  is  barely  detectable  (dashed 
line),  is  increasing  in  concentration  (double  lines),  is  becoming  de¬ 
pleted  (dotted  line),  or  is  in  relative  equilibrium  (single  line). 


Phase  I:  DO  Depletion  (DO  DEPLETE) 


68.  The  pH  range  over  which  the  phase  is  functional  is  6.3  to  7.5* 


69.  The  Eh  range  for  the  phase  is  +1*00  to  +600  mV. 
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TO.  The  DO  concentration  is  saturated  or  supersaturated  (at  given 
temperatures)  to  equal  approximately  2  mg/£  in  the  wa'oer  column. 

71.  Synopsis  of  phase.  Stratification  has  partially  isolated  the 
hypolimnion  from  diffusive  exchange  of  gases  with  surface  layers,  thus 
decreasing  reaeration.  The  topmost  layer  of  the  sediment  is  oxidized. 
Reduced  dissolved  organic  materials  and  ammonium  diffuse  from  the  under¬ 
lying  layers  of  sediment  to  the  top  layer;  however,  sufficient  oxygen 
remains  in  this  layer  to  oxidize  or  promote  microbial  oxidation  of  these 
substances  as  they  are  released.  Nitrate  diffuses  into  the  sediments 
from  the  overlaying  water  layer  and  is  denitrified  or  assimilated. 
However,  the  nitrification  process  in  the  hypolimnion  is  rapid  enough 

to  replenish  any  nitrate  lost  by  diffusion  or  flushing  of  the  system. 
Depletion  of  DO  exceeds  any  amount  that  might  be  brought  in  by  inflows. 
Inorganic  phosphate  is  continually  being  brought  in  by  inflows  and  is 
also  being  regenerated  from  organic  phosphate  within  the  system.  How¬ 
ever,  phosphate  is  also  lost  from  the  system  through  formation  and 
precipitation  of  ferric  phosphates.  Sulfate  may  be  accumulated  by  the 
system,  either  directly  through  inflows  bearing  high  levels  of  dissolved 
sulfate,  or  indirectly  via  oxidation  of  sulfides  released  by  decomposi¬ 
tion  of  organic  matter.  Carbon  dioxide  accumulation  results  from  res¬ 
piration  both  in  the  water  column  and  in  the  underlying  sediments. 

72.  Summary  of  processes.  The  components  and  pathways  for  Phase  I 
are  presented  in  Figure  11.  Arrows  in  Figure  11  are  explained  below: 


Arrow  No. 
1 

2 

3 

h 


_ Processes  Represented  by  Arrows _ 

Import  of  dissolved  inorganic  phosphate  with 
inflows . * 

Import  of  dissolved  inorganic  sulfate  with 
inflows . 

Import  of  dissolved  inorganic  nitrate  with 
inflows . 

Import  of  dissolved  inorganic  nitrite  with 
inflows . 

(Continued) 


*  Inflows  are  considered  major  movements  of  water  into  the 
hypolimnion  from  streams/rivers  that  flow  into  the  reser¬ 
voir  and  plunge  below  the  thermocline. 
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Arrow  No. 

5 

6 

7 

8 

9A 

9B 

10 

11 

12 

13 

14 

15 

15A 

15B 

15C 


15D 

15E 
1 6 


Processes  Represented  by  Arrows 


Import  of  dissolved  organic  and  particulate 
matter  (dissolved  organic  matter  with  inflows; 
particulate  matter  settling  from  metalimnion 
and  with  inflows ) . 

Import  of  ammonium  with  inflows. 

Export  of  dissolved  organics  and  particulate 
matter  with  outflows. 

Import  of  dissolved  carbon  dioxide  (predomi¬ 
nantly  as  HCO^)  with  inflows. 

Sedimentation  of  particulate  organic  matter  from 
this  water  layer  to  underlying  sediment  or 
water  layer. 

Release  of  dissolved  organic  matter  from  sedi¬ 
ment  (or  water  layer  beneath  layer  of  concern) 
to  this  water  layer  by  diffusion,  advection. 

Diffusion  of  ammonium  from  sediment  (or  under¬ 
lying  water  layer)  to  this  water  layer. 

Oxidation  of  ammonium  to  nitrite  (nitrification). 

Oxidation  of  nitrite  to  nitrate  (nitrification). 

Import  of  DO  with  inflows. 

Depletion  of  DO  by  diffusion  into  sediments  (or 
underlying  water)  layer. 


Oxygen  depletion  through  uptake  by  chemical 
oxygen  demand  and  biological  oxidation  pro¬ 
cesses,  which  may  be  split  into: 

0^  required  for  oxidation  of  NH^-^NOg. 

0g  required  for  oxidation  of  NO^NO^ . 

0g  required  for  oxidation  of  organically 
bound  phosphate  to  release  inorganic 
phosphate  (POjp. 

C>2  required  to  oxidize  organic  carbon  to 

->co2. 

02  required  to  oxidize  organic  sulfide  -♦■SO^. 

Decomposition  of  dissolved  organics  and  particu¬ 
late  matter  through  chemical  and  biological 
processes  to  yield  additional  dissolved  organ¬ 
ics  (retained  within  compartment),  and  carbon 
dioxide  and  dissolved  inorganic  phosphate, 
sulfate,  and  nitrate.  Process  may  be  split 
into: 


(Continued) 
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Arrow  No. 
l6A 

16b 

16c 

i6d» 

17 

18 

19 

20 


21 

22 

23 

2b 

25 


_ Processes  Represented  by  Arrows _ 

Oxidation  of  organically  1)01011  phosphate  to 
release  inorganic  phosphate  (POp. 

Oxidation  of  organic  sulfide/sulfate  to  re¬ 
lease  inorganic  sulfate  (SO^). 

Oxidation  of  organic  carbon  to  carbon  diox¬ 
ide  (co2). 

Oxidation  of  organic  nitrogen  to  inorganic 
nitrate  (NO^). 

Export  of  DO  with  outflows. 

Diffusion  of  carbon  dioxide  from  sediment  (or 
underlying  water)  layer  into  water  layer  of 
concern. 

Export  of  dissolved  carbon  dioxide  (as  HCO  ) 
with  outflows. 

Precipitation  of  dissolved  inorganic  phosphate 
through  formation  of  insoluble  ferric  and 
manganic  phosphates  and  by  sorption  to  other 
precipitating  materials. 

Export  of  dissolved  inorganic  phosphate  with 
outflows. 

Export  of  dissolved  inorganic  sulfate  with 
outflows. 

Diffusion  of  inorganic  nitrate  from  water  layer 
of  concern  into  sediment  (or  underlying  water) 
layer. 

Export  of  dissolved  inorganic  nitrate  with 
outflows. 

Diffusion  of  inorganic  nitrite  from  water  layer 
of  concern  into  sediment  (or  underlying  water 
layer ) . 


*  Differentiated  from  ammonification  under  aerobic 
conditions . 


73.  Algorithms  summarizing  the  various  processes  operative  on 
each  component  in  Figure  11  are  given  below: 
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Component 


d  DO/dt 


d  NH^/dt 


d  NO“/dt 


d  DOM*/dt 
d  COg/dt 


d  PO^/dt 


d  SO^/dt 


Factors  Influencing  Change  of 
_ Concentration  of  Component _ 

Dissolved 

Inflows  -  Outflows  -  DO  [NO^  Nitrate  Formation 

+  Dissolved  Organic  Matter  Decay  +  Particulate 

Organic  Decay  +  NH^  +  Nitrite  Formation  + 

Sediment  Decay  +  BOD  Demand  +  S  Sulfate 
Formation] . 

Inflows  +  Diffusion  from  Sediments  or  Underlying 

Anoxic  Water  Layer  +  Ammonification  -  NH^  [NH^ 

Nitrite  Formation  +  Biological  Uptake  +  Diffu¬ 
sion  into  Overlying  Water  Layer]. 

Inflows  +  NO^  [NH^  Nitrite  Formation]  +  NO^  [NO^ 

Nitrate  Reduction]  -  Inflows  -  NO^  [NO^ 

Nitrate  Formation  +  Diffusion  into  Sediments 
or  Underlying  Anoxic  Water  Layer  +  Diffusion 
into  Overlying  Water  Layer ] . 

Inflows  +  Sediment  Release  +  DOM  [Particulate 
Organic  Matter  Decay]  -  Outflows  -  DOM  [Dis¬ 
solved  Organic  Matter  Decay], 

Inflows  +  CO^  [Diffusion  Out  of  Sediments  or 

Underlying  Water  Layer  +  Dissolved  Organic 
Matter  Decay  +  Particulate  Organic  Matter 
Decay]  -  Outflows  -  CO^  [Carbonate  Formation 

+  Diffusion  into  Overlying  Water  Layer]. 
Inflows  +  PO^  [Dissolved  Organic  PO^  Decay  + 
Particulate  Organic  PO^  Decay  +  Cellular  Re¬ 
lease  of  PO^]  -  Outflows  -  POj^  [Formation  of 
Precipitates  and  Coprecipitates  +  Biological 
Uptake ] . 

Inflows  +  SO^  [Dissolved  Organic  Sulfate  Decay  + 
Particulate  Organic  Decay  +  S  Sulfate  Forma¬ 
tion]  -  Outflows  -  SO^  [Biological  Uptake], 

(Continued) 


*  DOM  =  dissolved  organic  matter. 


1*0 


[ 


Factors  Influencing  Change  of 

_ Concentration  of  Component _ 

Inflows  +  NO^  [NO^  Nitrate  Formation]  -  Outflows 
-  NO^  [Biological  Uptake  +  NO^  Reduction  + 
Diffusion  into  Sediments  or  Underlying  Anoxic 
Water  Layer]. 

Particulate 

Inflows  +  Settling  from  Overlying  Water  Layer  - 
Outflows  -  POM  [Particulate  Organic  Matter 
Decay  +  Settling  to  Underlying  Water  or  Sedi¬ 
ment  Layer]. 

Inflows  +  Settling  from  Overlying  Water  Layer  - 
Outflows  -  PIM  [Settling  to  Underlying  Water 
or  Sediment  Layer ] . 


*  POM  =  particulate  organic  matter  -  includes 
microorganisms . 

**  PIM  =  particulate  inorganic  matter. 


Component 
d  NO“/dt  = 


d  POM*/ dt  = 


d  PIM**/dt  = 


Phase  II:  DO  Exhaustion  (DO  EXHAUST) 

71*.  The  pH  range  over  which  the  phase  is  functional  is  6.3  to  7.5. 

75.  The  Eh  is  poised  at  approximately  +1*00  mV. 

76.  The  DO  concentration  is  greater  than  or  equal  to  0  mg/ i  and 
less  than  2  mg/£. 

77-  Synopsis  of  phase.  Dissolved  oxygen  concentration  in  the 
water  layer  has  dropped  well  below  the  2-mg/Jt  level.  Redox  is  poised 
at  =  +U00  mV  (at  pH  7-0)  by  the  N02~N0^  couple.  Sufficient  DO  is 
present  to  continue  oxidation  of  ammonium  to  nitrite  and  nitrate;  how¬ 
ever,  the  rate  is  much  diminished,  and  ammonium  is  present  in  readily 
detectable  quantities.  Moreover,  the  oxygen  depletion  rate  (arrows  lb, 
15,  and  17)  is  quite  rapid.  The  redox  will  not  drop  below  =  +  1*00  mV 
even  after  DO  has  reached  0  mg/ Si  as  long  as  a  sufficient  concentration 
of  NO^-NO^  remains.  As  in  Phase  I,  nitrate  can  still  enter  the  system 
with  inflows  and  is  still  produced  by  oxidation  of  ammonium  through 
nitrite.  However,  removal  of  nitrate  by  flushing  out  of  the  system 
and  by  diffusion  into  sediments  (or  underlying  anoxic  water  layer)  has 
become  appreciable  and  nitrate  can  no  longer  be  considered  as  being 


actively  accumulated.  Inorganic  sulfate  and  inorganic  phosphate  forma¬ 
tion  and  precipitation  as  iron  phosphate  remain  as  in  Phase  I. 

7 8.  Summary  of  processes.  Figure  12  presents  the  components  and 
pathways  of  importance  for  Phase  II.  These  are  identical  to  those  for 
Phase  I;  only  the  rates  have  changed.  See  Phase  I  (paragraph  72)  for 
description  of  processes  represented  by  arrows. 

79-  Algorithms  summarizing  the  various  processes  operative  on 
each  component  in  Figure  12  are  given  below: 


Component 


Factors  Influencing  Change  of 
Concentration  of  Components 


Dissolved 


d  DO/dt  =  Inflows  -  Outflows  -  DO  [NO^  Nitrate  Formation  + 

Dissolved  Organic  Matter  Decay  +  Particulate 

Organic  Matter  Decay  +  NH^  Nitrite  Formation  + 

Sediment  Decay  +  BOD  Demand  +  S  Sulfate 
Formation] . 

d  NH^/dt  =  Inflows  +  Diffusion  from  Sediments  or  Underlying 

Anoxic  Water  Layer  +  Ammonification  -  NH^  [NH* 

Nitrite  Formation  +  Biological  Uptake  +  Diffu¬ 
sion  into  Overlying  Water  Layer], 

d  NO^/dt  =  Inflows  +  NO^  [NH^  Nitrite  Formation]  +  NO^  [NO^ 

Nitrate  Reduction]  -  Outflows  -  NO^  [NO^ 

Nitrate  Formation  +  Diffusion  into  Sediments 
or  Underlying  Anoxic  Water  with  Ensuing  Denitri¬ 
fication  +  Diffusion  from  Overlying  Water 
Layer ] . 

d  NO^/dt  =  Inflows  +  NO^  [NOg  Nitrate  Formation]  -  Outflows 

-  NO^  [Biological  Uptake  +  NO^  Reduction  + 

Diffusion  into  Sediments  or  Underlying  Anoxic 
Water  Layer  with  Ensuing  Denitrification  + 
Diffusion  from  Overlying  Water  Layer]. 

d  DOM/dt  =  Inflows  +  Sediment  Release  +  DOM  [Particulate 
Organic  Matter  Decay]  -  Outflows  -  DOM  [Dis¬ 
solved  Organic  Matter  Decay]. 


(Continued) 


!6B 
16  A 


Figure  12.  Major  chemical  changes  and  pathways  associated  with  Phase  II  of  the 

DOSAP  algorithm 


Component 
d  COg/dt 


d  POj^/dt 


d  SOj^/dt 


Factors  Influencing  Change  of 
_ Concentration  of  Component _ 

Inflows  +  COg  [Diffusion  Out  of  Sediments  or 

Underlying  Water  Layer  +  Dissolved  Organic 
Matter  Decay  +  Particulate  Organic  Matter 
Decay]  -  CO^  [Outflows  +  Carbonate  Formation 

+  Diffusion  into  Overlying  Water  Layer]. 
Inflows  +  PO^  [Dissolved  Organic  PO^  Decay  + 
Particulate  Organic  PO^  Decay  +  Cellular  PO^ 
Release]  -  Outflows  -  POjj  [Formation  of  Pre¬ 
cipitates  and  Coprecipitates  +  Biological 
Uptake] . 

Inflows  +  SO^  [Dissolved  Organic  Decay  +  Partic¬ 
ulate  Organic  Decay  +  S  Sulfate  Formation]  - 
SO^  [Outflows  +  Biological  Uptake]. 

Particulate 


d  POM/dt  =  Inflows  +  Settling  from  Above  -  Outflows  -  POM 
[Particulate  Organic  Matter  Decay  +  Settling 
to  Below], 

d  PIM/dt  =  Inflows  +  Settling  from  Above  -  Outflows  -  PIM 
[Settling  to  Below]. 


Phase  III:  Ammonium 
Accumulation  (AMMON  ACCUM) 

80.  The  pH  range  over  which  the  phase  is  functional  is  6.6  to  T.k. 

81.  The  Eh  is  poised  at  approximately  +300  mV  (+h00  >  +200  mV). 

82.  The  DO  concentration  is  less  than  1.5  to  0  mg/l. 

83.  Synopsis  of  phase.  Dissolved  oxygen  has  been  nearly  or 
totally  depleted.  Nitrate  is  now  being  actively  removed,  and,  although 
some  nitrate  may  continue  to  be  introduced  in  inflows,  the  rate  of 
nitrate  input  (arrow  3)  is  exceeded  by  the  removal  in  outflows  (arrow 
21)  and  by  a  large  rate  of  diffusion  into  underlying  sediment  or  the 
anoxic  water  layer.  The  concentration  of  nitrate  is  sufficiently  low 
that  nitrate  reduction  (nitrate  respiration)  is  of  little  consequence. 
Ammonium  is  being  actively  accumulated  now  that  nitrification  is  un¬ 
available  as  a  sink.  Moreover,  as  oxygen  disappears,  ammonification 


of  organic  materials  in  the  water  column  (arrow  11 )  supplements  the 
diffusion  of  ammonium  from  the  sediments.  Ammonium  can  be  found  in 
waters  withdrawn  from  the  hypolimnion.  With  decreasing  redox  values  in 
the  total  absence  of  oxygen,  very  low  concentrations  of  reduced  manga¬ 
nese  begin  to  appear  in  the  water;  however,  at  redox  potentials  =  +300 
mV,  iron  and  manganese  tend  to  reoxidize  and  are  precipitated  out  of 
solution  (arrows  7  and  12). 

84.  Summary  of  processes.  Figure  13  gives  the  components  and 
pathways  of  importance  for  Phase  III.  Arrows  in  Figure  13  are  explained 
below : 

Arrow  No.  _ Processes  Represented  by  Arrows _ 

1  Import  of  dissolved  inorganic  phosphate  with 
inflows. 

2  Import  of  dissolved  inorganic  sulfate  with 
inflows . 


3  Import  of  dissolved  nitrate  with  inflows. 

4  Import  of  dissolved  and  particulate  organic 

matter  (dissolved  organic  matter  in  inflows; 
particulate  matter  settling  from  metalimnion 
and  with  inflows ) . 

5  Import  of  ammonium  with  inflows. 

6  Export  of  dissolved  organics  and  particulate 

matter  in  outflows. 

7  Import  of  dissolved  carbon  dioxide  (at  this  pH 

predominantly  HCO^)  with  inflows. 

8  Fermentative  decomposition  of  dissolved  and 

particulate  organic  matter  to  yield  organic 
acids  (which  for  now  are  retained  within  the 
compartment)  and  carbon  dioxide. 

9A  Sedimentation  of  particulate  organic  matter  from 

water  layer  of  concern  to  sediment  (or  water 
layer  beneath  it ) . 

9B  Release  of  dissolved  organic  matter  from  sedi¬ 

ment  (or  anoxic  water  layer  beneath  layer  of 
concern)  to  water  by  diffusion. 

10  Diffusion  of  ammonium  from  sediment  (or  anoxic 

water  layer  beneath  layer  of  concern)  to  water 


DOSAP  algorithm 


Arrow  No. 
11A 

11B 

12 

13 

lit 

15 

16 

IT 

18 

19 

20 

21 

22 

23 

2U 

25 


_ Processes  Represented  by  Arrows _ 

Ammonification  of  dissolved  and  particulate 
matter  containing  organic  nitrogen. 

Immobilization  of  ammonium. 

Export  of  ammonium  in  outflows. 

Diffusion  of  dissolved  CO^  (at  this  pH  range, 
mostly  HCO^)  from  sediments  or  underlying 

anoxic  water  layer  into  water  layer  of 
concern. 

Export  of  dissolved  CO in  outflows. 

Diffusion  of  dissolved  reduced  (mostly  chelated) 
iron  out  of  sediment  or  underlying  anoxic 
water  layer  into  water  layer  of  concern. 

Oxidation  and  precipitation  of  iron  from  water 
layer  of  concern  to  sediment  or  underlying 
anoxic  water  layer. 

Release  of  dissolved  inorganic  phosphate  from 
anaerobically  decomposing  dissolved  and 
particulate  matter  containing  organic 
phosphates. 

Export  of  dissolved  inorganic  phosphate  in 
outflows. 

Release  of  dissolved  inorganic  sulfate  from 
anaerobically  decomposing  dissolved  and 
particulate  organic  matter  containing  organic 
sulfates. 

Export  of  dissolved  inorganic  sulfate  in 
outflows. 

Export  of  dissolved  nitrate  in  outflowing 
waters. 

Diffusion  of  dissolved  nitrate  from  water  into 
underlying  sediment  or  anoxic  water  layer 
where  it  is  reduced. 

Diffusion  of  reduced  Mn++  from  sediments  or 
underlying  anoxic  water  layer  into  water  layer 
of  concern. 

Import  of  reduced  Mn++  with  inflowing  waters. 

Oxidation  and  precipitation  of  manganese  from 
water  layer  of  concern  to  sediment  or  under¬ 
lying  sediment  or  underlying  anoxic  water 
layer . 


1*7 


85.  Algorithms  summarizing  the  various  processes  operative  on 
each  component  in  Figure  13  are  given  below: 


Component  Factors  Influencing  Concentration  of  Components 

Dissolved 


d  O/dt  =  0.  Once  all  0^  has  been  exhausted,  the  solution 

to  the  equation:  Inflows  -  Outflows  -  DO  [N0^ 
Nitrate  Formation  +  Dissolved  Organic  Matter 
Decay  +  Particrlate  Organic  Matter  Decay  +  NH^ 
+  Nitrate  Formation  +  Sediment  Decay  +  BOD 

Demand  +  S  Sulfate  Formation]  =  0. 

d  NH^/dt  =  Inflows  +  Diffusion  from  Sediment  or  Underlying 
Anoxic  Water  Layer  +  Ammonification  -  Outflows 

-  NH^  [Biological  Uptake  +  Diffusion  into 

Overlying  Water  Layer]. 

d  NO^/dt  =  Inflows  +  NO^  [N0^  Reduction]  -  Outflows  -  NO^ 

[Diffusion  into  Sediments  or  Underlying  Anoxic 
Water  Layer  +  Diffusion  into  Overlying  Water 
Layer ] . 

d  NO^/dt  =  Inflows  -  Outflows  -  N0^  [Biological  Uptake  + 

N0^  Reduction  +  Diffusion  to  Underlying  Water 

or  Sediment  Layer]  +  Diffusion  from  Overlying 
Water  Layer. 

d  C02/dt  =  Inflows  +  CO^  [Diffusion  from  Sediments  or 

Underlying  Water  Layer  +  Dissolved  Organic 
Matter  Decay  +  Particulate  Organic  Matter 
Decay]  -  Outflows  -  CO^  [Carbonate  Formation 

+  Diffusion  to  Above  Water  Layer]. 

d  PO^/dt  =  Inflows  +  P0^  [Dissolved  Organic  PO^  Decay  + 

Particulate  Organic  P0^  Decay  +  Cellular 

Release  of  Pi]  -  Outflows  -  P0^  [Formation 

of  Precipitates  and  Coprecipitates  +  Bio¬ 
logical  Uptake]. 

d  S0^/dt  =  Inflows  +  SO^  [Dissolved  Organic  Sulfate  Decay 
+  Particulate  Organic  Sulfate  Decay]  -  Out¬ 
flows  -  S0^  [Biological  Uptake]. 

(Continued) 
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Component  Factors  Influencing  Concentration  of  Components 

d  DOM*/dt  =  Inflows  +  Sediment  Release  +  DOM  [Particulate 
Organic  Matter  Decay]  -  Outflows  -  DOM  [Dis¬ 
solved  Organic  Matter  Decay]. 

Particulate 

d  POM/dt  =  Inflows  +  Settling  from  Above  Water  Layer  - 
Outflows  -  POM  [Particulate  Organic  Matter 
Decay  +  Settling  to  Underlying  Water  or 
Sediment  Layer]. 

d  PIM/dt  =  Inflows  +  Settling  from  Above  Water  Layer  - 

Outflows  -  PIM  [Settling  to  Underlying  Water 
or  Sediment  Layer ] . 

*  Will  normally  be  present  in  waters  withdrawn  from 
hypolimnion. 

Phase  IV:  Manganese 
Accumulation  (MANG  ACCUM) 

86.  The  pH  range  over  which  the  phase  is  functional  is  6.6  to 

87.  The  Eh  is  +200  mV  +  +50  mV. 

88.  The  DO  concentration  is  0  mg/t. 

89.  Synopsis  of  phase.  Redox  has  dropped  to  =  +200  mV,  which  is 
the  changeover  point  from  oxidizing  to  reducing  conditions.  Both  dis¬ 
solved  ammonium  and  dissolved  reduced  manganese  are  being  activity  ac¬ 
cumulated  and  are  found  in  this  water  layer  if  it  is  being  withdrawn  as 
a  part  of  reservoir  releases.  Continued  withdrawal  of  this  water  will 
result  in  oxygen  sag  downstream  (classical  Streeter-Phelps )  due  to 
nitrification  of  ammonium  and  oxidation  of  reduced  manganese  to  yield 
manganic  hydroxides  and  oxides.  Reduced  iron  is  now  appearing  in  solu¬ 
tions  in  easily  detectable  concentrations,  although  most  of  this  will 
be  complexed  with  dissolved  organic  matter;  however,  at  a  redox  of 
+200  mV,  most  reduced  iron  is  chemically  reoxidized  to  yield  insoluble 
ferric  oxides  and  hydroxides  that  precipitate  back  out  of  the  water. 
Thus,  the  concentration  of  iron  in  solution  does  not  achieve  high 
levels.  Dissolved  nitrate  has  been  almost  totally  exhausted,  and  any 
residual  nitrate  is  totally  dependent  upon  the  amount  being  supplied  by 
external  sources,  i.e.  arrow  3.  Any  nitrate  entering  the  system  will 

L9 
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almost  immediately  oxidize  MN  to  MN  . 

90.  Summary  of  processes.  Figure  l4  depicts  the  major  components 
and  pathways  that  occur  in  Phase  IV.  Arrows  in  Figure  lU  are  explained 
below: 


Arrow  No. 
1 

2 

3 

h 


5 

6 

7 

8 


9A 

9B 

10 


11A 

11B 

12 

13 

lU 

15 

16 


_ Processes  Represented  by  Arrows _ 

Import  of  dissolved  inorganic  phosphate  with 
inflows . 

Import  of  dissolved  inorganic  sulfate  with 
inflows. 

Import  of  dissolved  nitrate  with  inflows. 

Import  of  dissolved  organic  and  particulate 
matter  (dissolved  organic  matter  in  inflows; 
particulate  matter  settling  from  metalimnion 
and  with  inflows ) . 

Import  of  ammonium  with  inflows. 

Export  of  dissolved  organics  and  particulate 
matter  with  outflows. 

Import  of  dissolved  carbon  dioxide  (as  HCO^). 

Fermentative  decomposition  of  both  dissolved  and 
particulate  organic  matter  to  yield  organic 
acids  (which  are  retained  within  this  compart¬ 
ment)  and  carbon  dioxide. 

Sedimentation  of  particulate  matter  from  water 
layer  of  concern  to  sediment. 

Release  of  dissolved  organic  matter  from  sedi¬ 
ment  to  overlying  water  by  diffusion. 

Diffusion  of  ammonium  from  sediment  or  anoxic 
water  layer  beneath  the  water  layer  of  concern 
to  overlying  water. 

Ammonif ication  of  dissolved  organic  and  particu¬ 
late  matter  containing  organic  nitrogen. 

Immobilization  of  ammonium. 

Export  of  ammonium  in  outflows. 

Diffusion  of  dissolved  CO^  (HCO^)  from  sediment. 

Export  of  dissolved  CO^  with  outflows. 

Diffusion  of  dissolved  reduced  (mostly  chelated) 
iron. 

Oxidation  and  pr^  ’’pitation  of  iron  from  water 
layer  of  concern  ^o  sediment. 

(Continued) 
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DOSAP  algorithm 


Arrow  No. 
17 


18 

19 


20 

21 

22 

23 

2k 

25 

26 


27 


_ Processes  Represented  by  Arrows _ 

Release  of  dissolved  inorganic  phosphate  from 
anaerobically  decomposing  dissolved  and  par¬ 
ticulate  matter  containing  organic  phosphates. 

Export  of  dissolved  inorganic  phosphate  with 
outflows. 

Release  of  dissolved  inorganic  sulfate  from 
anaerobically  decomposing  dissolved  and  par¬ 
ticulate  matter  containing  organic  sulfates. 

Export  of  dissolved  inorganic  sulfate  in  outflows 

Export  of  dissolved  nitrate  with  outflows. 

Diffusion  of  dissolved  nitrate  from  water  into 
underlying  sediment  where  it  is  denitrified. 

Diffusion  of  reduced  manganese  from  sediment 
into  water  layer  of  concern. 

Import  of  reduced  manganese  with  inflows. 

Export  of  dissolved  reduced  manganese  with 
outflows. 

Release  of  dissolved  reduced  manganese  from  anaer¬ 
obically  decomposing  dissolved  and  particulate 
matter  containing  organic  and  inorganic  Mn. 

Release  of  P0^  from  underlying  sediment  or 
anoxic  water  layer. 


91.  Algorithms  summarizing  the  various  processes  operative  on 
each  component  in  Figure  14  are  given  below: 


Factors  Influencing  Change  of 

Component  _ Concentration  of  Components _ 

Dissolved 

d  0/dt  =  0.  DO  concentration  =  0  mg/i. 
d  NH^/dt  =  Inflows  +  Diffusion  from  Sediment  or  Underlying 
Anoxic  Water  Layer  +  Ammonif ication  -  NH^ 
[Biological  Uptake  +  Diffusion  into  Overlying 
Water  Layer  +  Outflows*]. 

(Continued) 


*  Will  normally  be  present  in  waters  that  are  withdrawn 
from  the  hypolimnion. 
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Component 
d  NO“/dt 
d  NO~/dt 

d  DOM*/dt 
d  C02/dt 

d  PO^/dt 

d  SO^/dt 
d  Mn++/dt 


*  Normally 


Factors  Influencing  Change  of 
_ Concentration  of  Components _ 

0.  N02  concentration  =  0  mg/£. 

0.  If  >  0,  then  Inflows  -  Outflows  -  N0^  [Bio¬ 
logical  Uptake  +  N0^  Reduction  +  Diffusion 

into  Sediment  or  Underlying  Anoxic  Water 
Layer]  +  Diffusion  from  Overlying  Water 
Layer . 

Inflows  +  Sediment  Release  +  DOM  [Particulate 
Organic  Matter  Decay]  -  Outflows  -  DOM  [Dis¬ 
solved  Organic  Matter  Decay], 

Inflows  +  COg  [Diffusion  from  Sediments  or 

Underlying  Water  Layer  +  Dissolved  Organic 
Matter  Decay  +  Particulate  Organic  Matter 
Decay]  +  Carbonate  Formation  +  Diffusion  to 
Overlying  Water  Layer ] . 

Inflows  +  PO^  [Dissolved  Organic  PO^  Decay  + 
Particulate  Organic  P0^  Decay  +  Biological  P0^ 
Release]  -  Outflows  -  PO^  [Formation  of  Pre¬ 
cipitates  and  Coprecipitates  +  Biological  PO^ 
Uptake ] . 

Inflows  +  SO^  [Dissolved  Organic  Sulfate  Decay  + 

Particulate  Organic  Sulfate  Decay  -  Outflows  - 

S0^  [Biological  Uptake], 

Inflows  +  Mn++  [Diffusion  from  Underlying  Sedi¬ 
ments  or  Anoxic  Water  Layer  +  Dissolved  Or¬ 
ganic  Decay  +  Particulate  Organic  Decay  + 

Reduction  of  Mn  in  Suspended  Particulate 

Inorganic  Matter]  -  Outflows*  -  Mn++  [Bio¬ 
logical  Uptake  +  Diffusion  to  Overlying  Water 
Layer  +  Formation  of  Insoluble  Precipitates 

(C0^)  +  Oxidation  to  Mn+i+]. 


(Continued) 


inconsequential. 
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Factors  Influencing  Change  of 

Component  _ Concentration  of  Components _ 

d  Fe++/dt  =  Inflows  +  Fe++  [Diffusion  from  Underlying  Sedi¬ 
ment  or  Anoxic  Water  Layer  +  Dissolved  Organic 
Matter  Decay  +  Particulate  Organic  Matter  De¬ 
cay  +  Reduction  of  Fe+++  in  Suspended  Partic¬ 
ulate/Colloidal  Inorganic  Matter]  -  Outflows 

-  Fe++  [Biological  Uptake  +  Diffusion  to  Over- 
lying  Water  Layer  +  Formation  of  Insoluble 

Precipitates  (S  )  +  Oxidation  to  Fe+++]. 


Particulate 


d  POM/dt  =  Inflows  +  Settling  from  Overlying  Water  Layer  - 
Outflows  -  POM  [Particulate  Organic  Matter 
Decay  +  Settling  to  Underlying  Water  or  Sedi¬ 
ment  Layer]. 

d  PIM/dt  =  Inflows  +  Settling  from  Overlying  Water  Layer  - 
Outflows  -  PIM  [Settling  to  Underlying  Water 
or  Sediment  Layer ] . 


Phase  V :  Iron 
Accumulation  (FERR  ACCUM) 

92.  The  pH  range  over  which  the  phase  is  functional  is  6.6  to  7.1*. 

93.  The  Eh  is  +50  mV. 

9^.  The  DO  concentration  is  0  mg/ll. 

95-  Synopsis  of  phase.  Redox  has  dropped  to  +50  mV  and  reducing 
conditions  have  developed.  Iron  is  being  actively  reduced,  and  inor¬ 
ganic  phosphate  is  being  brought  into  solution  with  iron.  These  reduc¬ 
tion  products  are  found  in  the  outflows,  if  the  reservoir  has  hypolim- 
netic  withdrawal.  Ammonium,  manganese,  iron,  and  carbon  axe  being 
actively  accumulated  in  the  water,  although  the  rates  of  accumulation 
of  these  materials  are  somewhat  less  than  in  previous  phases;  these 
will  also  be  exported  in  the  outflows  from  a  system  that  has  hypolim- 
netic  withdrawal. 

96.  Summary  of  processes.  Figure  15  presents  the  components  and 
pathways  of  importance  for  Phase  V.  Arrows  in  Figure  15  are  explained 
below: 


Figure  15.  Major  chemical  changes  and  pathways  associated  with  Phase  V  of  the 

DOSAP  algorithm 


Arrow  No. 
1 

2 

3 

5 

6 
7 

8A 

8B 

9 

10A 

10B 

11 

12 

13 

14 


_ Processes  Represented  By  Arrows _ 

Import  of  dissolved  inorganic  phosphate  with 
inflows. 

Import  of  dissolved  inorganic  sulfate  with 
inflows. 

Import  of  dissolved  organics  and  particulate 
matter  (dissolved  organic  matter  in  inflows; 
particulate  matter  settling  from  metalimnion 
and  with  inflows ) . 

Import  of  ammonium  with  inflowing  waters. 

Export  of  dissolved  organics  and  particulate 
matter  with  outflows. 

Import  of  dissolved  carton  dioxide  (HCO  )  with 
inflows .  J 

Fermentative  decomposition  of  dissolved  organics 
and  particulate  organic  matter  to  yield 
organic  acids  (retained  within  dissolved 
organics  compartment)  and  dissolved  carton 
dioxide. 

Sedimentation  of  particulate  matter  from  the 
water  layer  of  concern  to  the  sed'-ment  layer 
beneath  it. 

Release  of  dissolved  organic  matter  from  sedi¬ 
ment  (water  layer)  beneath  the  water  layer  of 
concern  to  water  by  diffusion. 

Diffusion  of  ammonium  from  sediment  or  under¬ 
lying  anoxic  water  layer  into  water. 

Ammonification  of  dissolved  organics  and  por¬ 
tions  of  particulate  matter  containing  organic 
nitrogen. 

Immobilization  of  ammonium. 

Export  of  ammonium  in  outflows. 

Diffusion  of  dissolved  CO^  (HCO^)  from  sediments 
or  underlying  anoxic  water  layer  into  water. 

Export  of  dissolved  CO^  with  outflows. 

Release  of  dissolved  inorganic  phosphate  via 
decomposition  of  dissolved  organics  and  por¬ 
tions  of  particulate  organic  matter  containing 
phosphate. 


(Continued ) 


56 


Arrow  No. 
15 


_ Processes  Represented  by  Arrows _ 

Release  of  dissolved  inorganic  phosphate  from 
sediment  following  reductive  decomposition  of 
ferric  phosphates. 

1 6  Export  of  dissolved  inorganic  phosphate  with 

outflows . 

17  Release  of  reduced  iron  from  anaerobically  de¬ 

composing  dissolved  organics  and  particulate 
matter. 

18  Diffusion  of  dissolved  reduced  (mostly  chelated) 

iron  out  of  sediment  or  underlying  anoxic 
water  layer  into  water  layer  of  concern. 

19  Export  of  reduced  dissolved  iron  in  outflows. 

20  Release  of  dissolved  inorganic  sulfate  from 

anaerobically  decomposing  dissolved  organics 
and  particulate  matter. 

21  Export  of  dissolved  inorganic  sulfate  with 

outflows . 

22  Import  of  reduced  dissolved  manganese  with 

inflows . 

23  Release  of  reduced  dissolved  manganese  from 

anaerobically  decomposing  dissolved  organic 
and  particulate  matter. 

2k  Diffusion  of  reduced  dissolved  manganese  out  of 

sediment  or  underlying  anoxic  water  layer. 

25  Export  of  reduced  dissolved  manganese  with 

inflows . 


97.  Algorithms  summarizing  the  various  processes  operative  on 
each  component  in  Figure  15  are  given  below: 


Factors  Influencing  Change  of 

Component  _ Concentration  of  Components _ 

Dissolved 

d  0/dt  =  0.  DO  concentration  =  0  mg/l. 

d  NH^/dt  =  Inflows  +  Diffusion  from  Sediments  or  Underlying 
Anoxic  Water  Layer  +  Ammonification  -  Outflows 

-  NH^  [Biological  Uptake  +  Diffusion  into 

Overlying  Layer]. 

(Continued) 


Component 
d  NO"/dt 
d  NO~/dt 
d  DOM/dt 


d  C02/dt 


d  POjj/dt 


d  SO^/dt 
d  MN++/dt 


d  Fe++/dt 


I 


Factors  Influencing  Change  of 
_ Concentration  of  Components _ 

0.  N0~  concentration  =  0  mg/Jl. 

0.  NO”  concentration  =  0  rng/H. 

Inflows  +  Sediment  Release  (or  Underlying  Water 
Layer)  +  DOM  [Particulate  Organic  Matter 
Decay]  -  Outflows  -  DOM  [Dissolved  Organic 
Matter  Decay]. 

Inflows  +  C02  [Diffusion  from  Sediments  or 

Underlying  Water  Layer  +  Dissolved  Organic 
Matter  Decay  +  Particulate  Organic  Matter 
Decay]  -  Outflows  -  CO^  [Carbonate  Formation 

+  Diffusion  to  Overlying  Water  Layer]. 

Inflows  +  POj^  [Dissolved  Organic  P0^  Decay  + 

Particulate  Organic  P0^  Decay  +  Biological 

PO^  Release]  -  Outflows  -  POj^  [Formation  of 

Precipitates  +  Biological  PO^  Uptake]. 

Inflows  +  SO^  [Dissolved  Organic  Sulfate  Decay  + 

Particulate  Organic  Sulfate  Decay]  -  Outflows 

-  SO^  [Biological  Uptake]. 

Inflows  +  Mn++  [Diffusion  from  Underlying  Sedi¬ 
ments  or  Water  Layer  +  Dissolved  Organic 
Matter  Decay  +  Particulate  Organic  Matter 

Decay  +  Reduction  of  Mn  in  Suspended  Par- 
ticulate  Inorganic  Matter]  -  Outflows*  - 

Mn++  [Biological  Uptake  +  Diffusion  to  Over- 
lying  Water  Layer  +  Formation  of  Insoluble 

Precipitates**  +  Oxidation  to  Mn+1+]. 

Inflows  +  Fe++  [Diffusion  from  Underlying  Sedi¬ 
ment  or  Anoxic  Water  Layer  +  Dissolved  Organic 
Matter  Decay  +  Particulate  Organic  Matter 

Decay  +  Reduction  of  Fe+++  in  Suspended 
Particulate/Colloidal  Inorganic  Matter]  - 

Outflows*  Fe  [Biological  Uptake  +  Diffusion 
to  Overlying  Water  Layer  +  Formation  of  In¬ 
soluble  Precipitatest  +  Oxidation  to  Fe+++]. 


*  Present  in  outflows. 
**  Carbonates, 
t  Sulfides. 
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Factors  Influencing  Change  of 

Component  _ Concentration  of  Components _ 

d  S  /dt  =  Inflows  +  S  [Diffusion  from  Underlying  Sediment 
or  Anoxic  Water  Layer  +  Dissolved  Organic 
Matter  Decay*  +  Particulate  Organic  Matter 

Decay*  +  SO^  (Reduction  of  Sulfate)]  -  Out¬ 
flows  -  S"  [Biological  Uptake  +  Diffusion  to 
Overlying  Water  Layer  +  Formation  of  Insoluble 

Precipitates**  +  Oxidation  to  S°  or  SO^]. 
Particulate 

d  POM/dt  =  Inflows  +  Settling  from  Overlying  Water  Layer  - 
Outflows  -  POM  [Particulate  Organic  Matter 
Decay  +  Settling  to  Underlying  Water  or  Sedi¬ 
ment  Layer]. 

d  PIM/dt  =  Inflows  +  Settling  from  Overlying  Water  Layer  - 
Outflows  -  PIM  [Settling  to  Underlying  Water 
or  Sediment  Layer  +  Particulate  Inorganic 
Matter  Decomposition]. 

*  Sulfhydryl  groups  that  cannot  be  oxidized  to  S07  or  S°. 

**  Sulfides .  4 


Phase  VI:  Sulfide 
Accumulation  (SULFID  ACCUM) 

98.  The  pH  range  over  which  the  phase  is  functional  is  6.8  to  1.2. 

99.  The  Eh  is  approximately  -150  mV. 

100.  The  DO  concentration  is  0  mg/Jl. 

101.  Synopsis  of  phase.  Redox  has  descended  to  -150  mV.  Condi¬ 
tions  have  become  highly  reducing,  and  sulfide  release  from  sediments  or 
underlying  reduced  water  layers  has  begun.  Sulfate  has  begun  to  ac¬ 
tively  diffuse  into  sediment /underlying  anoxic  water  layer;  the  dis- 
similatory  sulfate  reduction  occurring  in  the  sediments  causes  the 
underlying  layers  to  serve  as  a  sink  for  sulfate.  Sulfate  continues 

to  move  into  the  system  via  inflowing  waters,  and  sulfate  is  also 
supplied  via  decomposition  of  organic  sulfates.  Ammonium,  reduced 
iron,  manganese,  and  phosphate  continue  to  be  released  from  the  sedi¬ 
ments  and  from  the  anaerobic  decomposition  of  organic  materials  bearing 
them.  Some  sulfide  has  begun  to  be  removed  by  formation  of  insoluble 


ferrous  sulfide.  However,  while  sulfide  is  actively  accumulating  in  the 
water  layer,  its  concentration  is,  as  yet,  too  low  for  significant  pre¬ 
cipitation  by  iron. 

102.  Summary  of  processes.  Figure  1 6  presents  the  components  and 
pathways  of  importance  for  Phase  VI.  Arrows  in  Figure  l6  are  explained 


below : 


I 


Arrow  No.  _ Processes  Represented  by  Arrows _ 

1  Import  of  dissolved  inorganic  phosphate  with 

inflows . 

2  Import  of  dissolved  inorganic  sulfate  with 

inflows . 

3  Import  of  dissolved  organics  and  particulate 

matter  (dissolved  organic  matter  in  inflows; 
particulate  matter  settling  from  metalimnion 
and  with  inflows ) . 

1+  Import  of  ammonium  with  inflows. 

5  Export  of  dissolved  organics  and  particulate 

matter  with  outflows. 

6  Import  of  dissolved  carbon  dioxide  (at  this  pH, 

mostly  HC0~ ) . 

T  Fermentative  decomposition  of  both  dissolved 

and  particulate  organic  matter  to  yield 
organic  acids  (retained  within  dissolved 
organics  compartment )  and  dissolved  carbon 
dioxide. 

8A  Sedimentation  of  particulate  matter  from  the 

water  layer  of  concern  to  the  sediment  (or 
anoxic  water  layer)  beneath  it. 

8b  Release  of  dissolved  organic  matter  from  sedi¬ 

ment  (or  anoxic  water  layer)  beneath  the 
water  of  concern  to  water  by  diffusion. 

9  Diffusion  of  ammonium  from  sediment  or  anoxic 

water  in  layer  beneath  the  layer  of  concern 
to  water  by  diffusion. 

10A  Aramonification  of  dissolved  organics  and  por¬ 

tions  of  particulate  matter  containing 
organic  nitrogen. 

10B  Immobilization  of  ammonium. 

11  Export  of  ammonium  in  outflowing  waters. 

(Continued) 
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Arrow  No. 
12 


13 
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15 


16 

17 


18 


19 

20 


21 

22 


23 

2k 

25 

26 


_ Processes  Represented  by  Arrows _ 

Diffusion  of  dissolved  CO^  (HCO^)  from  sedi¬ 
ments  or  underlying  anoxic  water  layer  into 
water. 

Export  of  dissolved  CO^  with  outflows. 

Release  of  dissolved  inorganic  phosphate  via 
decomposition  of  dissolved  organics  and  por¬ 
tions  of  particulate  matter  containing 
phosphate. 

Release  of  dissolved  inorganic  phosphate  from 
sediment  following  reductive  decomposition 
of  ferric  phosphates. 

Export  of  dissolved  inorganic  phosphate  with 
outflows . 

Release  of  reduced  iron  from  anaerobically  de¬ 
composing  dissolved  organics  and  particulate 
matter. 

Diffusion  of  dissolved  reduced  iron  out  of 
sediment  or  underlying  anoxic  water  layer 
into  water  layer  of  concern. 

Export  of  dissolved  reduced  iron  with  outflows. 

Release  of  sulfate  from  anaerobically  decom¬ 
posing  dissolved  organics  and  particulate 
matter . 

Export  of  dissolved  inorganic  sulfate  with 
outflows . 

Precipitation  of  reduced  dissolved  iron  to 
sediment  or  underlying  anoxic  water  layer  as 
a  consequence  of  the  formation  of  insoluble 
ferrous  sulfide. 

Diffusion  of  dissolved  inorganic  sulfate  into 
sediment  or  underlying  anoxic  water  layer. 

Diffusion  of  dissolved  inorganic  sulfide  out 
of  sediment  or  underlying  anoxic  water  layer. 

Export  of  dissolved  inorganic  sulfide  with 
outflows . 

Precipitation  of  dissolved  inorganic  sulfide 
to  sediment  or  underlying  anoxic  water  layer 
as  a  consequence  of  the  formation  of  insol¬ 
uble  ferrous  sulfide. 


(Continued) 
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27 

28 

29 


30 


_ Processes  Represented  by  Arrows _ 

Import  of  reduced  manganese  with  inflows. 

Release  of  reduced  manganese  from  underlying 
sediments  or  anoxic  water  layer. 

Release  of  reduced  manganese  from  anaerobically 
decomposing  dissolved  organics  and  particu¬ 
late  matter. 

Export  of  reduced  manganese  with  outflows. 


103.  Algorithms  summarizing  the  various  processes  operative  on 
each  component  in  Figure  l6  are  given  below: 


Factors  Influencing  Change  of 
Component  _ Concentration  of  Components 

Dissolved 


d  0/dt 
d  NH^/dt 


d  N0"/dt 
d  N0“/dt 
d  DOM/dt 


d  COg/dt 


d  P0^/dt 


0.  DO  concentration  =  0  mg/H. 

Inflows  +  Release  From  Sediment  or  Underlying 
Water  Layer  +  Ammonifieation  -  Outflows*  - 

NH^  [Biological  Uptake  +  Diffusion  into 

Overlying  Layer]. 

0.  N0~  concentration  =  0  mg/t. 

0.  N0^  concentration  =  0  mg/£. 

Inflows  +  Release  from  Sediment  or  Underlying 
Water  Layer  +  DOM  [Particulate  Organic  Matter 
Decay]  -  Outflows  -  DOM  [Dissolved  Organic 
Matter  Decay]. 

Inflows  +  CO ^  [Diffusion  from  Underlying  Sedi¬ 
ment  or  Water  Layer  +  Dissolved  Organic 
Matter  Decay  +  Particulate  Organic  Matter 
Decay]  -  Outflows  -  CO^  [Carbonate  Formation 

+  Diffusion  to  Overlying  Water  Layer]. 

Inflows  +  P0^  [Dissolved  Organic  PO^  Decay  + 
Particulate  Organic  P0^  Decay  +  Biological 
P0^  Release]  -  Outflows  -  POjj  [Formation  of 
Precipitates  +  Biological  P0^  Uptake]. 


*  Present  in  outflows. 
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Component 


Factors  Influencing  Change  of 
Concentration  of  Components 


d  SO^/dt  =  Inflows  +  SO^  [Dissolved  Organic  Sulfate  Decay 
+  Particulate  Organic  Sulfate  Decay]  -  Out¬ 
flows  -  SO^  [Biological  Uptake]. 

d  Mm++/dt  =  Inflows  +  Mn++  [Diffusion  from  Underlying  Sedi¬ 
ments  or  Water  Layer  +  Dissolved  Organic 
Matter  Decay  +  Particulate  Organic  Matter 

+4 

Decay  +  Reduction  of  Mn  in  Suspended  Par¬ 
ticulate  Organic  Matter]  -  Outflows*  -  Mn++ 
[Biological  Uptake  +  Diffusion  to  Overlying 
Water  +  Formation  of  Insoluble  Precipitates** 

+  Oxidation  to  Mn+^]. 

d  Fe++/dt  =  Inflows  +  Fe++  [Diffusion  from  Underlying  Sedi¬ 
ment  or  Anoxic  Water  Layer  +  Dissolved 
Organic  Matter  Decay  +  Particulate  Organic 

Matter  Decay  +  Reduction  of  Fe+++  in  Sus¬ 
pended  Particulate/Colloidal  Inorganic 

Matter]  -  Outflows*  -  Fe++  [Biological  Uptake 
+  Diffusion  to  Overlying  Water  Layer  +  Forma¬ 
tion  of  Insoluble  Precipitates**  +  Oxidation 

+++ 

to  Fe  ]. 

d  S  /dt  =  Inflows  +  S  [Diffusion  from  Underlying  Sedi¬ 
ment  or  Anoxic  Water  Layer  +  Dissolved 
Organic  Matter  Decayt  +  Particulate  Organic 

Matter  Decayt  +  SO^  Reduction]  -  Outflows*  - 

S  [Biological  Uptake  +  Diffusion  to  Over- 
lying  Water  Layer  +  Formation  of  Insoluble 

Precipitates**  +  Oxidation  to  S°/SO^]. 
Particulate 

d  POM/dt  =  Inflows  +  Settling  from  Overlying  Water  Layer  - 
Outflows  -  POM  [Particulate  Organic  Matter 
Decay  +  Settling  to  Underlying  Water  or 
Sediment  Layer ] . 

(Continued) 


*  Present  in  outflows. 

**  Sulfides. 

t  Sulfhydryl  groups  that  can  be  released  as  S_  or  SH 
cannot  be  oxidized  due  to  lack  of  0„. 
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Factors  Influencing  Change  of 

Component  _ Concentration  of  Components _ 

d  PIM/dt  =  Inflows  +  Settling  from  Overlying  Water  Layer  - 
Outflows  -  PIM  [Settling  to  Underlying  Water 
or  Sediment  Layer  +  Particulate  Inorganic 
Matter  Decomposition]. 


Phase  VII:  Methanogenesis 

10l*.  The  pH  range  over  which  the  phase  is  functional  is  6.6  to 

7.U. 

105.  The  Eh  is  less  than  or  equal  to  -175  mV. 

106.  The  DO  concentration  is  0  mg/£. 

107.  Synopsis  of  phase.  Redox  has  now  reached  -175  mV.  Ammonium, 
soluble  manganese  and  iron,  and  phosphate  continue  to  be  produced,  and 
these  materials  continue  to  appear  in  the  tailwaters  of  projects  with 
hypolimnetic  withdrawal.  In  those  impoundments  having  low  sulfide 
levels  in  their  sediments,  and  having  high  iron  levels,  the  rate  of 
sulfide  release  to  the  water  column  has  been  severely  diminished.  Con¬ 
sequently,  removal  of  sulfide  by  formation  of  insoluble  ferrous  sulfide 
and  by  export  from  the  system — either  in  outflows  or  by  the  activity  of 
sulfide-oxidizing  bacteria  at  the  metalimnetic-hypolimnetic  interface — 
results  in  a  net  depletion  of  sulfide  in  the  water  column.  Lack  of  sul¬ 
fate  as  an  electron  acceptor  has  also  caused  a  decrease  in  the  sulfide 
formation  rate.  The  decrease  in  activity  by  sulfate  reducers  has  re¬ 
moved  the  source  of  competition  for  organic  substrates  that  formerly 
existed  between  the  sulfate  reducers  and  the  methane-producing  bacteria 
and  methanogenesis  has  become  the  dominant  process.  As  a  result,  the 

Eh  has  shifted  from  -150  to  -175  mV  and  is  now  being  poised  by  methano¬ 
genesis.  Carbon  dioxide  is  being  actively  utilized  in  the  methane¬ 
forming  process. 

108.  Summary  of  processes.  Figure  17  depicts  the  components  and 
pathways  of  importance  for  Phase  VII.  Arrows  in  Figure  17  are  explained 
below: 


_ Processes  Represented  by  Arrows _ 

Import  of  dissolved  inorganic  phosphate  with 
inflows . 

Import  of  dissolved  and  particulate  organic 
matter  (dissolved  organic  matter  in  inflows; 
particulate  matter  settling  from  metalimnion 
and  with  inflows ) . 

Import  of  ammonium  with  inflows. 

Export  of  dissolved  organics  and  particulate 
matter  with  outflows. 

Import  of  dissolved  carbon  dioxide  (HCO  )  with 
inflows . 

Fermentative  decomposition  of  both  dissolved 
and  particulate  organic  matter  to  yield 
organic  acids  (retained  within  the  dissolved 
organics  compartment)  and  dissolved  carbon 
dioxide. 

Sedimentation  of  particulate  matter  from  the 
water  layer  of  concern  to  the  sediment  or 
anoxic  water  layer  beneath  it. 

Release  of  dissolved  organic  matter  from  the 
sediment  or  anoxic  water  layer  beneath  the 
layer  of  concern  into  the  water  in  the  layer 
of  concern  by  diffusion. 

Diffusion  of  ammonium  from  the  sediment  or 
anoxic  water  layer  beneath  the  layer  of  con¬ 
cern  into  the  water  in  the  layer  of  concern. 

Ammonium  release  by  ammonification  of  dissolved 
organics  and  portions  of  particulate  matter 
containing  organic  nitrogen. 

Ammobilization  of  ammonium. 

Export  of  ammonium  with  outflows  and  advection. 

Release  of  dissolved  CO^  (HCO^)  from  sediments 

or  underlying  anoxic  water  layer  into  the 
water  in  the  layer  of  concern  by  diffusion. 

Uptake  of  dissolved  CO^  (HCO_)  by 
methanogenesis. 

Export  of  dissolved  CO ^  (HCO.)  with  outflows 
and  advection. 


(Continued) 
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_ Processes  Represented  by  Arrows _ 

Release  of  dissolved  inorganic  phosphate  via 
fermentative  decomposition  of  dissolved 
organics  and  portions  of  particulate  matter 
containing  phosphate. 

Release  of  dissolved  inorganic  phosphate  from 
sediment  or  underlying  anoxic  water  layer  by 
diffusion  of  the  phosphate  after  reductive 
decomposition  of  ferric  phosphates. 

Export  of  dissolved  inorganic  phosphate  with 
outflows  and  advection. 

Release  of  dissolved  reduced  iron  from  the 
sediment  or  underlying  anoxic  water  layer  by 
diffusion  of  the  iron  after  reductive  decom¬ 
position  by  ferric  phosphates. 

Release  of  dissolved  reduced  iron  from  anaero¬ 
bically  decomposing  dissolved  organics  and 
particulate  matter. 

Precipitation  of  reduced  dissolved  iron  to 
sediment  or  underlying  anoxic  water  layer  as 
a  consequence  of  the  formation  of  insoluble 
ferrous  sulfide. 

Export  of  reduced  dissolved  iron  with  outflows 
and  advection. 

Production  of  methane  during  fermentation  of 
dissolved  organics  in  water  layer. 

Release  of  methane  to  the  water  layer  from 
sediment  or  underlying  anoxic  water  layer. 

Export  of  methane  from  system — primarily  as 
rising  bubbles  of  gas. 

Release  of  dissolved  reduced  manganese  from 
underlying  sediment  or  anoxic  water  layer 
by  diffusion  into  the  layer  of  concern. 

Import  of  reduced  manganese  with  inflows  and 
advection. 

Release  of  reduced  manganese  from  anaerobically 
decomposing  dissolved  organics  and  particu¬ 
late  matter. 

Export  of  reduced  manganese  with  outflows  and 
advection. 


(Continued) 
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_ Processes  Represented  by  Arrows _ 

Precipitation  of  dissolved  sulfide  to  sediment 
or  underlying  anoxic  water  layer  as  a  conse¬ 
quence  of  the  formation  of  insoluble  ferrous 
sulfide. 

Release  of  dissolved  sulfide  from  sediment  or 
underlying  anoxic  water  layer  to  water  layer 
of  concern  by  diffusion. 

Export  of  dissolved  sulfide  with  outflows  and 
advection. 


109.  Algorithms  summarizing  the  various  processes  operative  on 
each  component  in  Figure  17  are  given  below: 


Component 


d  0/dt 
d  NH^/dt 


d  N0~/ dt 
d  N0”/dt 
d  DOM/dt 


d  C02/dt 


d  POjj/dt 


Factors  Influencing  Change  of 
_ Concentration  of  Components _ 

Dissolved 

0.  DO  concentration  =  0  mg/i. 

Inflow  +  Release  from  Sediment  or  Underlying 
Water  Layer  +  Ammonification  -  Outflows  - 

NH^  [Biological  Uptake  +  Diffusion  into 

Overlying  Layer]. 

0.  N02  concentration  =  0  mg/i.. 

0.  N0^  concentration  =  0  rng/J l. 

Inflows  +  Release  from  Sediment  or  Underlying 
Water  Layer  +  DOM  [Particulate  Organic  Matter 
Decay]  -  Outflows  -  DOM  [Dissolved  Organic 
Matter  Decay], 

Inflows  +  C02  [Diffusion  from  Underlying  Sedi¬ 
ment  or  Water  Layer  +  Dissolved  Organic 
Matter  Decay  +  Particulate  Organic  Matter 
Decay]  -  Outflows  -  CO^  [Carbonate  Formation 

+  Diffusion  to  Overlying  Water  Layer]. 

Inflows  +  PO^  [Dissolved  Organic  P0^  Decay  + 
Particulate  Organic  P0^  Decay  +  Biological 
P0^  Release]  -  Outflows  -  P0j^  [Formation  of 
Precipitate  +  Biological  P0^  Uptake], 

(Continued) 
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Component 
d  SO^/dt  = 


Factors  Influencing  Change  of 
Concentration  of  Components 


Inflows  +  SO^  [Dissolved  Organic  Sulfate  Decay 
+  Particulate  Organic  Sulfate  Decay]  -  Out- 


++ 


flows  -  SO^  [Biological  Uptake]. 

++ 


Mn'/dt  =  Inflow  +  Mn''  [Diffusion  from  Underlying  Sedi¬ 
ments  or  Water  Layer  +  Dissolved  Organic 
Matter  Decay  +  Particulate  Organic  Matter 

+1* 

Decay  +  Reduction  of  Mn  in  Suspended  Par¬ 
ticulate  Organic  Matter]  -  Outflows  -  Mn++ 
[Biological  Uptake  +  Diffusion  to  Overlying 
Water  +  Formation  of  Insoluble  Precipitates 

+  Oxidation  to  Mn+^ ] . 

d  Fe++/dt  =  Inflows  +  Fe++  [Diffusion  from  Underlying  Sedi¬ 
ment  or  Anoxic  Water  Layer  +  Dissolved  Or¬ 
ganic  Matter  Decay  +  Particulate  Organic 

Matter  Decay  +  Reduction  of  Fe+++  in  Sus¬ 
pended  Particulate/Colloidal  Inorganic 

Matter]  -  Fe++  [Biological  Uptake  +  Diffusion 
or  Overlying  Water  Layer  +  Formation  of  In¬ 
soluble  Precipitates  +  Oxidation  to  Fe+++], 

d  S  /dt  =  Inflows  +  S  [Diffusion  for  Underlying  Sediment 
or  Anoxic  Water  Layer  +  Dissolved  Organic 
Matter  Decay  +  Particulate  Organic  Matter 

Decay  +  SO^  Reduction]  -  Outflows  -  S  [Bio¬ 
logical  Uptake  +  Diffusion  to  Overlying  Water 
Layer  +  Formation  of  Insoluble  Precipitate  + 


Oxidation  to 


SO"]. 


d  CH^/dt  =  CH^  [Diffusion  from  Underlying  Sediment  or  Anoxic 
Water  Layer]  +  CH^  [Methanogenesis  (Dissolved 
Organic  Matter  Decay  +  Particulate  Organic 


Matter  Decay)]  -  Outflows  - 
Diffusion  to  Overlying  Water  Layers]. 


CH^  [Oxidation  + 


Particulate 

d  POM/dt  =  Inflows  +  Settling  from  Overlying  Water  Layer  - 
Outflows  -  POM  [Particulate  Organic  Matter 
Decay  +  Settling  to  Underlying  Water  or 
Sediment  Layer]. 

(Continued) 


Component 


Factors  Influencing  Change  of 
Concentration  of  Components 


d  PIM/dt  =  Inflows  +  Settling  from  Overlying  Water  Layer  - 
Outflows  -  PIM  [Settling  to  Underlying  Water 
or  Sediment  Layer  +  Particulate  Inorganic 
Matter  Decomposition]. 


PART  V:  EVALUATION  OF  RELEVANT  RATE  DATA 


110.  Little  information  exists  in  the  literature  on  rate  data 
for  anaerobic  processes  in  reservoirs  included  in  DOSAP.  By  contrast, 
there  is  considerable  information  on  these  processes  in  marine  and 
freshwater  ecosystems  other  than  reservoirs.  Rather  than  furnishing  a 
compendium  of  nonreservoir  related  data,  a  tabulation  is  presented 
herein  of  ranges  of  values  that  are  specific  for  reservoirs  (Table  l); 
values  for  nonreservcir  systems  are  included  only  when  there  are  nc 
data  for  particular  processes  in  reservoirs  (Table  l).  For  convenience 
in  making  comparisons  with  data  from  the  authors  own  studies  and  in 
order  to  establish  a  common  mode  for  data  assembly,  all  of  this  informa¬ 
tion  is  expressed  in  terms  of  milligrams  of  the  component  being  released 

-2  -1 

or  depleted  per  square  metre  of  surface  per  day  (mg*m  *day  ). 

111.  While  extensive  use  has  been  made  of  the  literature  on 
flooded  soils  as  a  source  of  concepts  for  the  DOSAP  algorithm,  this  area 
has  not  been  used  as  a  source  of  rate  information.  There  are  several 
reasons  for  this.  Many  of  the  flooded  soils  experiments  that  have  been 
done  were  conducted  at  temperatures  in  the  25°  to  35°C  range.  Since 
temperatures  in  reservoir  hypolimnions  typically  are  in  the  5°  to  15°C 
range,  the  rates  obtained  from  these  studies  tend  to  be  higher  than 
anticipated  for  reservoir  situations.  Moreover,  many  flooded  soils 
experiments  tend  to  use  large  amendments  of  the  particular  substrate 
being  studied,  and  this  may  also  yield  above-normal  rates.  Also,  the 
previous  history  of  many  test  soils  with  respect  to  flooding  is  not 
described  in  the  literature  in  which  the  experiments  are  presented. 

Thus,  the  applicability  of  data  obtained  using  these  soils  is  somewhat 
dubious;  the  authors  cannot  be  certain  whether  the  rate  data  obtained 

in  such  studies  would  be  typical  of  newly  flooded  soils,  soils  located 
in  the  flood  pool  region  of  a  reservoir,  or  reservoir  sediments. 

Finally,  experiments  on  flooded  soils  often  use  a  very  thin  layer  of 
water  over  the  soil  surface.  This  practice  results  in  a  larger  potential 
capacity  for  diffusive  exchange  of  oxygen  throughout  the  entire  layer 
of  water  over  the  soil  surface  than  would  be  expected  in  the  typical 


reservoir  situation.  The  larger  potential  capacity  would  tend  to 
produce  a  higher  concentration  of  oxidized  chemical  species  in  the  water 
over  the  flooded  soil  than  would  be  anticipated  for  all  but  the  shal¬ 
lowest,  most  well -mixed  reservoirs. 

112.  It  is  apparent  from  the  data  presented  in  Table  1  that  there 
is  some  information  available  on  almost  all  of  the  anaerobic  release 
processes  of  DOSAP.  Inflow  and  outflow  rates  are  necessarily  site 
specific  for  a  given  reservoir  and  will  not  be  considered  herein.  How¬ 
ever,  because  of  the  restricted  amount  of  data  available  for  reservoirs, 
there  are  a  number  of  limitations  on  what  may  be  done  with  the  data. 

For  example,  the  flux  rates  presented  in  Table  1  are  a  useful  indication 
of  the  range  of  concentrations  of  anaerobic  products  that  may  be  expected 
in  reservoirs;  however,  no  indications  are  given  of  the  rates  of  releases 
to  be  expected,  and  rates  are  mandatory  for  DOSAP  to  be  used  as  the 
basis  for  a  mathematical  subroutine  simulating  anaerobic  processes. 
Moreover,  the  data  in  Table  1  are  also  limited  relative  to  the  number 
of  environmental  situations  to  which  they  may  be  applied.  For  example, 
the  anoxic  reservoir  hypolimnion  data  presented  in  Table  1  were  taken 
from  one  environment  and  the  temperature  range  in  this  case  was  10°  to 
lU°C;  the  coastal  marine  sediment  data  were  taken  from  a  i*°C  environ¬ 
ment;  and  the  denitrification  data  obtained  from  the  sediment  of  the 
small  eutrophic  lake  were  obtained  at  a  third  temperature.  The  question 
that  immediately  arises  is  what  would  the  values  in  Table  1  be  if  they 
were  all  obtained  from  the  same  system  at  the  same  temperature.  To 
accomplish  this,  however,  studies  must  be  conducted  under  highly  con¬ 
trolled  conditions. 

113.  The  data  presented  in  Table  2*  illustrate  the  use  of  sedi¬ 
ment  water  reaction  chambers  to  provide  data  for  DOSAP.  While  the  data 
in  Table  2  are  expressed  as  fluxes  (net  rate  and  direction  of  release), 
they  could  also  have  been  given  as  first-order  rate  coefficients  or 
other  desired  means  of  expression.  The  central  point  here  is  that  all 

*  Table  2  presents  the  results  of  studies  conducted  by  the  authors  on 
sediment  water  reaction  in  DeGray  Reservoir,  Arkansas. 


the  rate  data  necessary  for  DOSAP  can  be  collected  using  this  procedure. 
To  provide  a  complete  set  of  rate  coefficients  for  DOSAP  under  a  wide 
range  of  geographic  conditions,  the  authors  are  studying  the  anaerobic 
release  of  these  same  components  from  sediments  taken  from  various  CE 
reservoir  projects  in  different  sections  of  the  continental  United 
States.  Rate  coefficients  obtained  from  these  studies  will  be  summarized 
in  an  Engineering  Technical  Letter  {ETL)  on  data  preparation  and  coeffi¬ 
cient  selection  (September  1982).  This  ETL  will  also  give  detailed  in¬ 
structions  on  the  methods  to  be  used  in  collecting  the  data  (or  for 
making  estimates,  when  necessary)  required  for  DOSAP.  For  illustrative 
purposes.  Table  3  presents  a  compendium  of  the  methods  used  to  obtain 
data  required  for  use  in  the  DOSAP  model  and  is  intended  as  preliminary 
guidance  prior  to  publication  of  the  ETL. 


PART  VI:  CONCLUSIONS 


llU.  The  occurrence  of  DO  depletion  and  anaerobic  processes  in 
the  hypolimnion  of  CE  reservoirs  can  adversely  impact  water  quality  and 
directly  affect  CE  Civil  Works  activities. 

115.  A  review  of  the  literature  and  discussions  with  authorities 
on  anaerobic  processes  confirmed  the  feasibility  of  formulating  a 
realistic  conceptual  model  of  anaerobic  processes  that,  along  with 
appropriate  process  rate  data,  could  serve  as  a  basis  for  a  subroutine 
for  numerical  reservoir  water  quality  models.  A  numerical  model  with 
the  ability  to  simulate  deoxygenation  and  anaerobic  processes  is 
required  to  provide  the  CE  with  the  ability  to  address  problems  associ¬ 
ated  with  low  DO  and  anaerobic  conditions  in  reservoirs. 

116.  Some  of  the  data  required  to  use  DOSAP  will  be  site  specific 
and  pertinent  field  data  must  be  collected  for  each  reservoir  applica¬ 
tion  to  provide  a  description  of  initial  conditions  and  to  update  input 
variables.  However,  for  developing  and  verifing  DOSAP  as  an  anaerobic 
subroutine  for  widespread  field  application,  a  broad  base  of  data 
obtained  from  a  variety  of  different  reservoirs  is  required.  Sediment 
water  reaction  chambers  offer  the  most  viable  means  for  obtaining  such 
information.  These  required  studies  are  being  accomplished  as  part  of 
the  Corps'  EWQOS  Program. 

117.  Application  of  the  DOSAP  subroutine  in  a  numberical  water 
quality  model  to  various  CE  reservoir  projects  using  data  obtained  from 
the  projects  offers  a  viable  means  for  examining  the  effectiveness  of 
various  management  strategies  upon  problems  associated  with  low  DO  and 
anaerobic  conditions.  Under  proper  conditions  and  through  the  use  of 
data  obtained  from  neighboring  bodies  of  water,  a  water  quality  model 
equipped  with  the  DOSAP  subroutine  should  also  be  applicable  to  preim¬ 
poundment  investigations  that  involve  a  concern  for  the  probability  of 
a  new  project  to  develop  anoxic  conditions. 
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_ Component _  _ Method* _  _ References _ 

Dissolved  oxygen  Biochemical  oxygen  demand  probe  with  associated  DO  meter;  Standard  Methods  (1976) 

calibrated  by  lodometric  Method  'for  DO  measurement 

Oxidation-reduction  potential  Platinum  electrode  and  standard  pH  calomel  reference  electrode  Platinum  electrode  preparation 
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